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Abstract 

Several  shell  deformation  models  are  developed  for  use  in  crashworthiness  analysis  of  rotationally 
symmetric  structures.  These  models  use  analytical  techniques  to  predict  the  crushing  force 
versus  axial  crush  distance  characteristics  of  both  a  rigid-plastic,  hemispherical  shell  and  an 
elastic,  cylindrical  shell  loaded  axially  by  a  rigid  flat,  plate.  Additional  methods  are  proposed  to 
determine  the  effects  of  cutout  sections  and  internal  stiffening  members  on  the  crushing  force 
capacity  of  the  shells.  These  proposed  methods  are  applied  to  determine  the  energy  absorption 
capability  of  the  composite/metal  nose  structure  of  a  mini  submarine  subjected  to  a  head-on 
impact  with  a  flat  rigid  wall.  The  nose  structure  is  composed  of  a  stiffened  composite  shell  that 
covers  and  is  attached  to  the  metal,  forward,  hemispherical  portion  of  the  pressure  hull.  The 
complex  structure  is  simplified  to  a  rotationally  symmetric  shell  model  and  modified  to  account 
for  the  effects  of  stiffening  elements  and  cutout  sections.  Laminated  plate  theory,  a  progressive 
composite  failure  method,  and  the  models  developed  are  then  used  to  determine  the  structure's 
energy  absorption  capacity.  The  fairing  structure  is  determined  to  be  capable  of  absorbing  the 
kinetic  energy  associated  with  an  impact  from  an  initial  vehicle  speed  of  7.5  knots.  Further,  the 
entire  nose  structure  (the  fairing  and  forward  pressure  hull)  is  predicted  to  be  able  to  absorb 
the  kinetic  energy  from  a  26.4  knot  vehicle  impact. 
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Thesis  Supervisor:  James  H.  Williams,  Jr. 

Title:  SEPTE  Professor  of  Engineering 
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Chapter  1 


Introduction 


The  Advanced  Seal  Delivery  System  (ASDS)  (Figure  (1-1))  is  a  highly  sophisticated,  mini¬ 
submarine  19.8  m  long  that  displaces  60  metric  tons  submerged.  It  is  designed  to  transport  a 
team  of  Navy  SEALs  from  a  host  submarine  to  a  tactical  point  of  insertion  and  back.  While 
the  transported  SEALs  are  a  formidable  offense  weapon,  the  ASDS  itself  has  no  offensive  or 
defensive  capabilities  and  must  rely  on  stealth  and  its  ability  to  ’’hide”  in  the  water  column 
to  ensure  the  safe  transport  of  its  personnel.  As  such,  the  ASDS  will  be  operating  in  close 
proximity  to  the  ocean  bottom  in  order  to  mask  its  underwater  signature.  This  near  bottom 
operating  strategy  presents  some  problems,  however.  The  most  pressing  of  these  is  the  increased 
risk  of  impact  with  submerged  objects  or  structures.  For  the  safety  of  the  embarked  crew,  it  is 
paramount  that  a  determination  be  made  as  to  the  crashworthiness  of  the  ASDS  and  whether 
operational  restrictions  need  to  be  imposed.  In  effect,  a  safety  envelope  needs  to  be  developed 
which  takes  into  account  the  depth,  bottom  type  and  impact  angle  to  give  a  maximum  safe 
operating  speed.  A  natural  starting  point  for  this  analysis  is  to  look  at  how  a  front  end  impact 
will  affect  the  forward  structure  of  the  vessel.  This  analysis,  however,  presents  no  small  task. 
The  complexity  of  this  problem  quickly  becomes  apparent  when  the  structure  types  and  problem 
variables  are  viewed: 

•  The  structure  of  the  nose  section  consists  of  various  composite  materials  and  steel  in-, 
eluding:  glass  woven  cloth,  carbon  cloth,  kevlar.  vinylester,  rubber,  syntactic  foam,  and 
HY-80  steel. 
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Figure  1-1:  External  Isometric  View  of  the  ASDS 

•  The  variables  of  the  problem  include:  depth,  impact  angle,  speed,  and  bottom  type. 

Additionally,  the  interaction  of  the  structural  material  and  modes  of  failure  during  defor¬ 
mation  and  fracture  are  not  well  understood  or  predictable  to  a  high  degree  of  accuracy.  In 
light  of  these  complexities,  using  model  testing  to  develop  some  empirical  means  of  relating 
the  structural  behavior  during  a  collision  to  key  parameters  such  as  impact  angle,  velocity,  or 
crush  distance  would  be  the  optimal  analysis  tool.  However,  due  to  the  cost  associated  with 
such  tests,  this  method  is  not  feasible.  Moving  to  the  other  end  of  the  spectrum,  a  computer 
simulated  ’’crash*'  test  of  the  ASDS  could,  in  theory,  be  used  to  ascertain  the  desired  safety 
envelope  at  a  substantially  cheaper  price  and  in  an  accelerated  time  frame.  While  modeling 
the  full  ASDS  via  computer  simulation  seems  plausible  on  the  surface,  a  little  digging  reveals 
that  even  the  most  sophisticated  computer  tools  have  not  evolved  far  enough  to  accurately 
predict  impact  responses  for  such  complex,  composite/metal  structures.  With  this  in  mind,  it 
is  apparent  that  neither  model  testing  nor  computer  simulation  alone  can  solve  this  problem. 
A  middle  road  between  these  methods  may  be  the  required  path.  In  effect,  a  combination  of 
structural  simplification  and  computer  and  analytical  modeling  will  be  used  to  complete  the 
analysis.  With  all  of  this  said,  the  overall  objective  of  the  project  will  be  to  determine  a  safe 
impact  speed  for  the  vehicle  such  that  during  an  impact  the  pressure  hull  will  not  rupture  and 
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the  crew  will  be  able  to  survive  the  collision. 

As  a  first  step  in  solving  the  problem,  this  thesis  will  analyze  the  energy  absorbing  potential 
of  the  composite  nose  fairing  and  forward  pressure  hull  hemi-head  structures  during  a  head-on 
impact  with  a  rigid  wall  structure.  While  this  impact  scenario  is  idealized  and  likely  never 
to  occur  during  actual  operations,  it  will  serve  as  a  worst  case  impact  and  best  case  energy 
absorption  scenario.  This  will  be  a  worst  case  impact  from  the  stand  point  that  all  kinetic  energy 
will  be  absorbed  by  the  structure  in  the  form  of  deformation.  No  energy  will  be  absorbed  by 
the  impacted  structure  or  be  maintained  as  vehicle  kinetic  energy.  This  is  a  best  case  energy 
absorption  scenario  in  that  the  impact  will  allow  the  full  utilization  of  the  bow  structure  to 
absorb  energy  before  the  pressure  hull  is  impacted.  The  reasoning  behind  looking  at  such  a 
case  first,  is  to  ascertain  the  maximum  energy  absorption  potential  of  the  composite  structure. 

This  analysis  will  progress  in  three  phases: 

•  Deformation  model  development. 

•  Structural  modeling  of  the  ASDS  nose  structure. 

•  Energy  absorption  analysis. 

The  models  developed  include  deformation  relationships  for  a  rigid-plastic,  hemispherical 
and  an  elastic,  cylindrical  shell  loaded  axially  by  a  rigid  flat  plate.  A  method  is  also  included  to 
approximate  the  reduced  load  carrying  capacity  caused  by  cutouts  in  the  shell  structure.  Next, 
the  composite  shell  of  the  nose  fairing  will  be  modeled  using  simplified,  rotationally  symmetric 
shells  (i.e.- 'cylinders,  frustum,  and  hemispheres).  Internal  stiffening  members  are  incorporated 
into  the  model  by  “smearing”  the  effective  area  of  each  stiffener  across  portions  of  the  inner 
surface  of  the  shell  to  increase  the  effective  thickness  of  the  structure.  Cutouts  in  the  shell  for 
components,  such  as  thruster  doors,  video  equipment,  and  lights  are  also  modeled.  Finally,  the 
deformation  relationships  developed  will  be  applied  to  the  corresponding  sections  of  the  model 
to  produce  a  force  verse  displacement  relationship  of  the  structure  as  it  is  crushed  axially  by  a 
rigid  flat  plate.  The  energy  absorption  capability  of  the  entire  structure  will  then  be  determined 
and  translated  into  the  vehicle  speed  needed  to  produce  such  an  impact. 
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Chapter  2 

Development  and  Discussion  of 
Analysis  Tools 

2.1  Crushing  Response  of  a  Hemispherical  Shell  Between  Rigid 
Plates 

2.1.1  Prior  Research 

Much  attention  has  been  devoted  to  the  study  of  the  energy  absorbing  characteristics  of  thin- 
walled,  metal,  dome  structures  [l]-[7] .  The  problem  of  large  deformations  of  rigid-plastic,  spher¬ 
ical  shells  being  compressed  between  rigid  plates  was  first  studied  by  Updike  [1],  A  simplified 
yield  condition  and  geometry  was  used  to  obtain  a  sequence  of  limit  loads,  thus  producing  a 
force-displacement  relationship  valid  from  a  displacement  of  a  few  shell  thicknesses  to  about 
1  /  10th  the  shell  radius.  Wierzbicki  [4]  used  similar  geometric  and  deformation  assumptions 
to  derive  a  force-deflection  relationship  based  on  the  optimization  of  bending  and  membrane 
energies  valid  theoretically  until  the  plate  deflection  reaches  half  the  radius  of  the  shell.  Com¬ 
parison  of  [1]  and  [4]  show  that  the  force-displacement,  results  predicted  by  each  differ  only 
slightly.  Experimental  testing  of  hemispherical  shells  was  reported  by  Kitching  [5]  and  Kin- 
head  [5] .  Kitching  et  al  studied  the  quasi-static  loading  of  hemispherical  shells  and  observed 
experimentally  three  different  stages  of  deformation  for  shells  with  radius  to  thickness  ratios 
ranging  from  36  to  460.  Later,  Kinkead  performed  experimental  testing  on  shells  with  radius 


11 


y* 


v 


Figure  2-1:  Geometry  of  Hemispherical  Deformation  Model. 

to  thickness  values  ranging  from  8  to  32  and  identified  two  vice  three  deformation  stages.  The 
loading  characteristic  observed  by  Kinhead  were  compared  to  values  predicted  by  the  methods 
of  [1]  and  [4]  and  large  differences  were  found.  Gupta  [7]  developed  another  prediction  method 
using  a  rigid  plastic  analysis  which  relies  on  two  calibration  constants  that  must  be  determined 
from  curve  fits  of  experimental  data.  The  method  proposed  by  Gupta  showed  good  correlation 
with  experimental  results.  This  analysis  develops  a  means  of  predicting  the  load-compression 
relationship  of  hemispherical  shells  using  a  rigid  plastic  analysis  technique  with  an  emphasis 
on  creating  a  simple,  analytical  model  whose  accuracy  is  improved  over  that  proposed  in  [1] 
and  [4].  The  model  developed  will  utilize  the  same  deformation  model  used  by  [1]  and  [4],  but 
calculate  the  membrane  energy  using  a  simple  strain  averaging  method.  Additionally,  the  re¬ 
sulting  relationship  will  be  capable  of  analyzing  material  with  differing  compressive  and  tensile 
stresses. 

2.1.2  Deformation  Process 
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The  spherical  shell  of  radius  R  will  be  loaded  from  the  top  by  a  moving  rigid  plate  with 
the  base  of  the  shell  pinned  to  a  fixed  rigid  plate.  As  the  shell  is  loaded  axially,  a  spherical 
dimple  forms  which  is  an  inversion  of  the  original  shell  curvature  as  shown  in  Figure  (2-1).  In 
this  deformed  state,  surface  AB  remains  rigid  and  undeformed.  At  point  B,  a  circumferential 
plastic  hinge  increases  the  curvature  of  the  shell,  followed  by  rigid  body  translation  between  B 
and  D  and  a  reversal  of  curvature  by  the  plastic  hinge  at  D.  The  shell  is  subsequently  unloaded 
and  rigid  body  translation  results  on  surface  DE.  As  the  shell  deforms,  the  hinge  circles  at  B 
and  D  move  progressively  outward.  Deformation  continues  in  this  manner  until  the  inverted 
portion  of  the  shell  at  point  E  contacts  the  fixed  rigid  plate. 

2.1.3  Deformation  Model 
Geometry 

This  analysis  begins  with  a  look  at  the  various  geometrical  relationships  that  can  be  attained 
from  the  deformation  model  in  Figure  (2-1).  First,  the  radius  to  the  center  of  the  toroidal 
surface  (a)  and  the  displacement  of  the  moving  rigid  plate  (S)  can  be  expressed  in  terms  of  the 
shell  radius  ( R )  and  deformation  angle  (<f>)  as 

a  =  i2sin0,  and  (2.1a) 

6  =  R(l—cos4>)  (2-lb) 

By  differentiating  Equations  (2.1a)  and  (2.1b),  the  radial  and  axial  velocities  of  the  toroidal 
section  and  rigid  plate  become 

da 

dt 
dS 

dt 

By  approximating  cos  <fr(t)  using  the  first  two  terms  of  its  Taylor  expansion  in  Equation  (2.1b) 
and  rearranging,  a  simple  equation  relating  the  deformation  angle  to  the  rigid  plate  displacement 
can  be  deduced. 

(2.4) 


R( cosd>)^y-  ,  and 
dt 


=  R  (sin  <f) ) 


do 

dt 


(2.2) 

(2.3) 
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Turning  now  to  the  toroidal  surface  BCD,  the  average  plastic  hinge  length  (L)  and  toroidal 
surface  area  (A)  can  be  calculated  by  the  Equations  (2.5a)  and  (2.5b). 


L  =  27ra,  and 
A  =  Anratp 


(2.5a) 

(2.5b) 


Additionally,  the  radial  and  rotational  velocities  of  the  plastic  hinge  at  point  B  can  be  taken 


de  =  v^Rdt 

dt  r  r  dt 

Finally,  in  order  to  simplify  subsequent  calculation,  surface  BCD  will  be  approximated  as  a 
parabolic  surface  of  the  form 

y==^7^c:x2+r^~cos^  (2-8) 

1  r  sm  cp 

with  the  x  and  y  axes  oriented  as  shown  in  Figure  (2.2.3). 

Derivation  of  Membrane  Stresses  and  Fully  Plastic  Bending  Moment 

Next,  the  membrane  force  ( N )  and  fully  plastic  bending  moment  (M)  relationships  for 
a  material  having  differing  compressive  ( ac )  and  tensile  (at)  flow  stresses  will  be  developed. 
As  shown  in  Figure  (2-2),  when  ac  ^  at  the  neutral  axis  will  shift  away  from  the  geometrical 
center  of  the  cross  section  by  an  amount  rj.  By  using  the  fact  that  the  neutral  axis  of  a  specimen 
subjected  to  a  bending  moment  is  located  where  the  membrane  force  (N)  is  equal  to  zero,  we 
are  able  to  calculate  a  general  expression  for  r)  in  terms  of  shell  thickness  (t),  ac,  and  cr(  as 


shown  below 


acdz  =  0 


where  both  at  and  ac  are  taken  to  be  positive..  Upon  integrating  and  solving  for  tj.  Equation 
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Figure  2-2:  Membrane  Stress  and  Fully  Plastic  Bending  Moment  Representation  for  Material 
having  Differing  Compressive  and  Tensile  Flow  Stresses. 


(2.9)  becomes 


1  jrc-<*t 

V  —  tt - 

2  Gc  +  (Tf 


With  rj  calculated,  general  expressions  for  N  and  M  can  be  derived  as  follows: 


(2.10) 


f3+»?  r2-V 

M  —  /  atzdz  +  /  aczdz 

Jo  Jo 


Note  that  the  regions  of  membrane  deformation  are  compressive  in  nature  while  the  bending 
forces  are  both  tensile  and  compressive.  After  integrating,  substituting  Equation  (2.10)  for  77, 
and  setting  a  =  2  ,  the  total  membrane  force  and  fully  plastic  bending  moments  become 


N 

M 


act 


<n 

<Tc  +  <Tt 


(Tct2 


4 


a 


(2.11) 

(2.12) 


In  order  verify  the  validity  of  Equations  (2.11)  and  (2.12),  crc  can  be  set  equal  to  at  which 
results  in  a  symmetric  stress  profile  throughout  the  specimen.  In  such  a  case  Equations  (2.11) 
and  (2.12)  reduce  to  the  standard  expressions  for  the  fully  plastic  membrane  force  and  bending 
moment. 
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Energy  Balance 


The  principle  of  virtual  velocities  will  be  used  to  relate  the  rate  at  which  energy  is  being 
introduced  to  the  system  to  the  rate  at  which  it  is  being  absorbed.  Equation  (2.13)  states  that 
the  rate  at  which  energy  is  applied  to  the  shell  must  exactly  match  the  rate  at  which  it  is  being 
absorbed  by  the  bending  and  membrane  effects  in  the  shell. 


(i)  (2)  (3) 

dd  _dE\j  dEm 
dt  dt  +  dt 


(2.13) 


The  remainder  of  this  section  will  further  develop  each  of  the  terms  in  Equation  (2.13).  Starting 
first  with  the  left  hand  side  of  Equation  (2.13),  term  (1)  can  be  rewritten  in  a  different  form 
by  using  Equation  (2.3)  as  follows 

PM~PRm*t  (2'14) 

The  second  term  of  Equation  (2.13)  represents  the  rate  at  which  energy  is  absorbed  by  the  shell 
due  to  bending  processes.  Using  the  deformation  model  shown  in  Figure  (2-1),  it  is  plainly  seen 
that  bending  energy  will  be  absorbed  at  two  discrete  locations  -  plastic  hinges  B  and  D.  This 
component  can  be  quantified  by  the  following  expression 

it  -  E  *§*  <2-I5a) 

i=B,D 

Equation  (2.15a)  can  be  further  simplified  by  making  the  following  assertions 

•  the  magnitudes  of  the  bending  moment  and  rate  of  angular  rotation  at  hinges  B  and  D 
are  equal. 

•  the  total  hinge  length  ^  L,  can  be  rewritten  as  2 L  (where  L  is  the  average  hinge 

i=B,D 

length. 
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Figure  2-3:  Deformation  of  Surface  BCD. 


Using  these  simplifications  and  Equations  (2.1a),  (2.5a),  and  (2.7),  Equation  (2.15a)  can  be 
expressed  in  the  following  form: 


dEb 

dt 


=  4ttR2  sin 


(2.16) 


The  final  term  in  Equation  (2.13)  accounts  for  the  membrane  energy  dissipated  on  surface  BCD 
and  can  be  written  as: 

d-¥  =  NATt  <2-17> 
where  ^  is  the  instantaneous  rate  of  strain  experienced  by  surface  BCD.  As  a  simplification, 
the  average  strain  rate,  will  be  used  in  place  of  the  instantaneous  strain  rate.  Figure 
(2-3)  shows  a  detailed  view  of  the  deforming  surface  BCD.  The  average  strain  of  this  surface 
will  be  calculated  by  dividing  the  surface  into  two  subsurfaces,  calculating  the  average  strain 
on  each,  and  then  averaging  the  strains  from  each  region  to  get  the  final  average  strain.  This 
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statement  can  be  expressed  mathematically  as 


Eav  — 


£Iav  +  Slla 


(2.18) 


where  ejav  and  ejjav  are  the  average  strains  experienced  as  surfaces  BC0  (region  1)  and  C0D0 
(region  2)  deform  to  BC and  CD  respectively.  Looking  first  at  region  I,  the  strain  can  be  written 
in  terms  of  initial  and  final  curvatures  of  the  BC. 


1  ( dy\ 2  1  2 

£l  =  £l}  ~  £l°  =  2  \dx )  ‘2  + 


(2.19a) 


In  general,  the  average  strain  in  region  I  can  be  formulated  as: 

rr  sin  <p 


£i« 


i  r 

r  sin  <j>  J0 


£jdx 


(2.20) 


Upon  substitution  of  Equation  (2.19a),  integration  and  simplification,  Equation  (2.20)  reduces 
to 

(2.21) 


=  Ia? 


£Iav  =  3^' 

The  strain  in  region  II  requires  an  additional  step  to  calculate.  While  the  strain  in  region  I 
is  purely  compressive  for  the  entire  deformation,  the  strain  in  region  II  is  initially  compressive 
for  the  first  part  of  deformation  and  then  tensile  for  the  remainder  of  the  deformation.  So, 
the  average  strain  will  be  calculated  for  the  compressive  phase  as  surface  C0D0  is  deformed  to 
an  intermediate  flat  surface  C,  A  and  also  for  the  tensile  phase  as  the  shell  further  deforms  to 
surface  CD.  With  this,  the  average  strain  in  region  II  can  be  calculated  in  the  same  manner 
as  in  region  I.  The  details  of  the  calculation  follow. 


£IIav 

^ ^avCOmpressiVe  ^ ^avtensitc 

^  rr  sirup 

r  / 

\ 

/ 

£Ilav 

r  sin  4>  Jo 

j  Uuf 

_  Y  J  compressive 

—  Ctt  I  +  [  £■ 

°compressive  J  V  J 

£IIav 

j  rr  sin  <p 

{b2-°)+ 

(-S' 

f— vv 

dx 

r  sin  o  J0 

\dx J  J 

£IIav 

-  b2 

dx 


(2.22) 
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With  the  average  strain  calculated  in  both  regions,  we  can  now  determine  the  average  strain 
of  the  entire  surface  BCD  by  substituting  Equations  (2.21)  and  (2.22)  into  Equation  (2.18). 
Subsequently,  the  average  strain  rate  can  be  calculated  by  simple  differentiation  of  the  average 
strain  as  follows. 


\<?  (2-23) 

4  <2-24> 

Now  a  final  expression  for  the  rate  at  which  energy  is  absorbed  by  membrane  forces  can  be 
realized  by  substituting  Equations  (2.5b),  (2.1a)  and  (2.24)  into  Equation  (2.17). 

^11  =  Anr NR  sin  4>62^-  (2.25) 

at  at 


deav 

dt 


New  expressions  have  now  been  formulated  for  each  term  in  Equation  (2.13).  By  substituting 
the  new  expressions  (Equations  (2.14),  (2.16),  and  (2.25))  into  Equation  (2.13)  and  solving  for 
the  instantaneous  crushing  force  (P),  the  following  result  is  obtained 


P  =  4tt 


^RM  +  Nr^_^ 


(2.26) 


Further,  the  mean  crushing  force  (P)  can  be  calculated  and  minimized  with  respect  to  the 
toroidal  radius  (r)  to  obtain  an  optimal  toroidal  radius. 


(2.27) 


(2.28) 


By  substituting  the  optimal  toroidal  radius  back  into  Equation  (2.27)  and  using  the  fact  that 
at  the  optimal  radius  the  membrane  and  bending  contributions  are  equal,  the  mean  crushing 
force  can  be  expressed  more  simply  as 
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(2.29) 


Popt  =  4  <pR 


Equation  (2.29)  can  be  rewritten  in  terms  of  plate  displacement  and  shell  thickness  by  using 
Equations  (2.4),  (2.11),  and  (2.12). 


Finally,  the  optimal  instantaneous  crushing  force  becomes 


(2.30) 


Pi Opt  =  4nV2acaH2S2 


opt 


(2.31) 


2.2  Derivation  of  Cylindrical  Tube  Equation 

2.2.1  Prior  Research 

The  deformation  and  energy  absorption  characteristics  of  thin-walled,  composite,  cylindrical 
shells  have  been  studies  for  many  years  by  many  researchers.  [8]-[l5].  Farley  [10]-[13]  and  Hull 
[14]  have  performed  extensive  testing  on  composite  shell  samples  composed  of  various  mate¬ 
rials  using  different  lay-up  conditions.  Both  quasi-static  and  dynamic  loading  situations  were 
investigated.  These  studies  did  much  to  further  the  knowledge  base  in  the  areas  of  determining 
crushing  modes  and  their  controlling  mechanisms.  They  also  showed  that  there  was  no  explicit 
set  of  rules  that  could  be  generalized  for  all  composites  types.  Material  lay-up  pattern,  crushing 
speed,  fiber  volume  fraction,  and  radius  to  laminate  thickness  ratio  all  play  key  roles  in  deter¬ 
mining  the  deformation  characteristics  of  the  shell.  From  these  studies,  complicated  numerical 
models  have  been  established  for  use  in  predicting  the  failure  of  very  simple  composite  struc¬ 
tures.  The  degree  of  success  realized  by  these  models  has  varied  in  even  simple  situations.  The 
vast  majority  of  the  work  described  has  been  directed  toward  understanding  the  deformation 
of  composites  on  the  micro-  and  macroscopic  scales.  Hanefi  and  Wierzbicki  [15]  took  a  more 
global  approach  to  a  similar  deformation  problem.  They  developed  a  simple  model  that  was 
used  to  predict  the  energy  absorption  characteristics  of  a  combined  metallic/composite  tube 
structure  in  axial  compression.  The  model  coupled  a  rigid  plastic  analysis  of  the  metal  tube 
with  a  compressive  perfectly  plastic  and  tensile  elastic  analysis  of  the  composite  material.  The 
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Figure  2-4:  Deformation  Model  for  Cylindrical  Thin-walled  Structure. 

predicted  results  matched  well  with  experimental  data.  In  this  analysis,  a  global  failure  ap¬ 
proach  will  also  be  taken.  A  deformation  model  will  be  developed  using  elastic  principles  to 
determine  the  force  deflection  and  energy  deflection  characteristics  of  a  thin-walled,  cylindrical 
shell  loaded  axially.  Bending  and  in-plane  stiffnesses  will  be  step-wise  varied  to  account  for 
stiffness  variations  as  plys  fail  within  the  structure.  Finally,  a  simple  result  will  be  established 
for  a  material  exhibiting  constant  and  identical  flexural  and  inplane  stiffness  properties. 

2.2.2  Deformation  Process 

Figure  2-4  shows  a  representation  of  two  buckling  lobes  of  the  cylinder  during  deformation 
where  6  denotes  the  axial  crush  distance,  a  is  the  radius  of  the  cylinder,  and  2 H  is  the  half 
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length  of  the  buckling  lobe  generated.  As  the  cylinder  is  crushed  axially  by  a  distance  8,  the 
cylinder  will  form  an  axisymmetric,  elastic  buckling  pattern  with  a  wave  length  equal  to  2 H. 
Deformation  takes  place  both  in  the  form  of  bending  and  circumferential  strain.  The  material 
will  behave  elastically  until  its  elastic  limits  are  reached  with  subsequent  deformation  occurring 
with  the  material  in  a  progressively  degraded  state.  The  material  degradation  takes  place 
through  a  step-wise  reduction  of  the  material’s  flexural  and  circumferential  elastic  moduli  and 
material  thickness.  The  process  continues  progressively  through  all  buckling  lobes  formed. 

2.2.3  Deformation  Model 

The  principle  of  virtual  work  (Equation  (2.32))  will  be  used  as  the  starting  relationship  from 
which  to  derive  an  expression  relating  the  crushing  force  ( P )  to  the  axial  crush  distance  (8). 
Equation  (2.32)  states  that  the  energy  applied  to  the  cylinder  through  an  external  force  (1) 
must  be  exactly  equal  to  the  energy  absorbed  by  the  structure  in  the  form  of  circumferential 
(2)  and  bending  (3)  strains. 


(1) 

f  (2) 

(3) 

PdS  = 

J  ( (Jgdeg  + 

axdex)dV 

(2.32) 

ae  = 

Eeg  and  crx 

=  Ee x 

(2.33) 

Pd6  = 

/  ( Eegdeg  +  Eexdex)dV 

Jv 

(2.34) 

Because  the  deformation  will  be  elastic,  Equation  (2.32)  can  be  simplified  using  the  elastic 
relationships  between  stress  and  strain  (Equation  (2.33))  to  obtain  Equation  (2.34).  Figure  (2- 
5)  shows  a  detailed  view  of  the  deformation  of  a  single  buckling  lobe  with  parameters  pertinent 
to  the  analysis  labeled.  Using  this  figure,  expressions  can  be  derived  in  terms  of  the  variables 
shown  for  the  incremental  volume  (dV),  strains  (ex  and  eg),  and  incremental  strains  (dsx  and 
deg).  First,  inextensibility  will  be  assumed  for  the  shell  as  it  buckles,  making  dV  constant  over 
the  deformation  process  and  having  the  following  value: 

dV  =  AiratH  (2.35) 
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Figure  2-5:  Detailed  View  of  the  Deformed  Cylinder  Section. 
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Next,  the  circumferential  strain  can  be  expressed  as  the  ratio  of  the  outward  deflection  of  the 
shell  ( w )  to  the  original  shell  radius.  The  actual  deflection  is  measured  from  the  original  surface 
of  the  shell  to  the  deformed  surface  and  varies  from  0  at  the  being  of  the  lobe  to  a  maximum 
value  of  w0  at  the  midpoint,  and  then  back  to  0  at  the  end  of  the  lobe.  In  order  to  simplify 
future  calculations,  a  linear  deflection  profile  ( w(x  )  =  jjw0)  is  used  to  approximate  the  actual 
deflection  profile  as  shown  in  Figure  (2-5).  Using  this  linear  relationship,  eg  and  deg  can  be 
expressed  as 


eg  = 

W  X 

o  =  ~aHWo 
/ 

(2.36) 

deg  = 

hdw° 

(2.37) 

The  bending  strain  component  can  be  determined  by  using  the  curvature  of  the  shell  section 
(k)  and  the  position  within  the  thickness  of  the  material  (z)  as  is  shown  in  Equation  (2.38). 
This  variation  of  the  bending  strain  throughout  the  shell  thickness  will  be  assumed  linear  with 
maximums  occurring  at  |  and  — |  and  no  strain  occurring  at  the  mid-plane.  Further,  plane 
sections  will  remain  plane  as  shell  curvature  increases. 


(2.38) 


However,  the  radius  of  curvature  ( R )  is  related  to  the  half  lobe  length  and  the  maximum 
circumferential  deflection  by  Equation  (2.39) 


R2  -  (R-  w0 )2  =  H2  -  w2 

H 2 

R  =  —  (2.39) 

2w0  v  ’ 

Substitution  of  Equation  (2.39)  into  Equation  (2.38)  yields  the  final  expressions  for  ex  and  dsx 
(Equations  (2.40)  and  (2.41)). 

(2.40) 

(2.41) 
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Now,  Equations  (2.35)-(2.37)  and  (2.40)  can  be  substituted  back  into  Equation  (2.34)  to  obtain 
the  following: 


Pd6  = 


w0dw0  (2.42) 


If  the  flexural  and  membrane  Young’s  moduli  can  be  assumed  equal  and  constant  throughout 
the  deformation  process,  Equation  (2.42)  can  be  further  simplified  to  Equation  (2.43  ). 


Pd6  =  -? tE 

O 


tH  at3] 

T  +  i? 3 


WndWn 


(2.43) 


Previous  work  [18]  has  proven  that  the  elastic,  buckling  lobe  length  of  a  cylinder  can  be  related 
to  the  shell  radius  and  thickness  by 


2  H  =  1.72  Vat 
H  »  fat 

Using  Equation  (2.44),  Equation  (??)  can  be  reduced  to 


(2.44) 


(2.45) 


A  final  simplification  can  be  made  to  Equation  (2.45)  by  relating  6  to  wQ.  Using  Figure  (2-5) 
the  following  geometric  relationships  can  be  generated. 


wn 


sm(a)  =  -jj 
cos  (a)  =  1~2 Jj 


(2.46) 
.  (2-47) 
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By  squaring  and  adding  Equations  (2.46)  and  (2.47),  the  relationships  are  combined  in  the 
following  form 


1 


ru>ol2 

r  s  i 

y + 

i 

d 

i _ 

1  2 


(2.48) 


After  differentiating,  rearranging  and  substituting  Equation  (2.44),  the  desired  relationship  is 
achieved 


wodw0 


dS 

2 


This  expression  can  now  be  substituted  back  into  Equation  (2.45)  to  obtain  the  final  expression 
relating  the  crushing  force  applied  to  the  axial  deflection  of  the  cylinder. 


(2.49) 


As  previously  states,  Equation  (2.49)  applies  to  materials  with  equal  and  constant  flexural  and 
circumferential  Young's  moduli.  If  these  conditions  are  not  met,  the  more  general  expression 
(Equation  (2.42))  must  be  used. 


2.2.4  Correction  for  Open  Sections 

Open  sections  or  cutouts  in  a  shell  will  effectively  reduce  the  load  required  to  deform  the  shell 
in  the  local  region  of  the  cutout.  These  cutouts  could  even  change  the  local  deformation  mode 
of  the  structure  or  change  the  order  in  which  the  overall  structure  deforms.  For  simplicity,  this 
analysis  will  include  only  the  reduction  in  loading  capacity  for  several  reasons.  The  relationships 
previously  developed.  Equations  (2.26)  and  (2.49),  are  valid  only  for  the  failure  modes  assumed; 
so,  a  change  in  failure  modes  (i.e.-  from  axisymmetric  shell  yielding  to  lobar  buckling  for  ex¬ 
ample)  will  not  be  captured  here  in.  While  this  seems  like  a  fairly  unreasonable  assumption, 
the  reader  is  reminded  that  the  objective  of  this  analysis  is  to  determine  the  maximum  energy 
absorption  capability  of  the  structure.  The  development  of  Equations  (2.31)  and  (2.49)  assume 
failure  modes  which  produce  this  maximum.  So,  a  shift  to  a  different  mode  would  produce  a 
lower  energy  absorption.  With  the  assumption  that  the  failure  mode  remains  unchanged  at  the 
cutouts,  it  does  not  make  sense  to  consider  the  order  in  which  parts  of  the  structure  collapse 
because  the  same  amount  of  energy  will  be  absorbed  regardless  of  the  collapse  sequence.  So,  the 
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only  effect  that  will  be  considered  is  the  reduction  of  strength  caused  by  the  loss  of  material. 
This  effect  will  be  approximated  by  multiplying  the  applied  load  of  an  intact  shell  by  a  reduc¬ 
tion  factor  in  the  local  vicinity  of  the  cutout.  This  reduction  factor  will  be  proportional  to  the 
remaining  shell  area  divided  by  the  material  present  in  an  intact  shell  of  the  same  dimensions 
(Equation  (2.50)) 


F(6)corrected  — 
e(6)  = 


F (6)  uncorrectedE  (6)  where 
(^complete  shell  (<5)  A-removed  material  (^)) 


■^complete  shell  (^) 


(2.50) 


While  this  area  ratio  is  part  of  the  reduction  effect,  it  does  not  capture  all  aspects.  Other 
factors,  such  as  cutout  geometry,  affect  the  load  carrying  capacity  of  the  structure  around 
the  cutout.  While  these  other  effects  are  recognizable,  a  theoretical  derivation  of  their  effect  is 
impractical.  This  effect  must  be  determined  experimentally.  The  experimental  correction  factor 
would  be  used  to  match  experimental  test  data  with  theoretical  predictions  and  is  expressed 
as  £  in  Equation  (2.50).  Since  no  experimental  calibration  was  done  during  the  course  of  this 
thesis,  £  will  be  set  equal  to  1  until  future  evaluation  determines  its  correct  value. 


Hemispherical  Dome 

For  the  hemispherical  dome,  the  cutout  sections  will  be  assumed  to  be  circular  with  diameter 
D  as  shown  in  Figure  (2-6).  From  Figure  (2-6)  it  is  apparent  that  the  correction  factor  can  be 
stated  in  terms  of  the  average  plastic  hinge  radius  (a(<5))  and  angle  <p(6)  as  follows: 


e(6) 

s(6) 


(Average  Hinge  Diameter- Length  of  Open  Section)  \  & 

1  S  hemisphere 


2ira(6) 


Average  Hinge  Diameter 
j>(6)a(8) 


2tt  a(S) 


^  £ hemisphere  ^  ^  £ hemisphere  (2-51) 


The  task  now  becomes  to  express  (f)(6)  in  terms  of  known  quantities.  Figures  (2-7)  and  (2-8) 
show  more  detailed  views  of  the  side  and  top  of  the  shell  respectively  with  pertinent  dimensions 
labeled  in  each  figure. 
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Figure  2-7:  Geometrical  Relationship  of  Parameters  Used  in  the  Derivation  of  the  Correction 
Factor  for  Open  Sections  in  a  Dome  Structure  (side  view). 

From  these  figures,  the  following  geometric  relationships  can  be  established: 


Al 

A2 

A3 

A4 

A5 


r0  +  D  cos  9  —  a(6 ) 
Al 
cos  9 


D-A2 

j-A3  =  A2-% 


(2.52) 


(f>(6)  can  now  be  expressed  in  terms  of  the  initial  radius  at  which  the  cutout  is  encountered  (r0), 
cutout  diameter  (D),  cutout  inclination  angle  (9),  and  average  plastic  hinge  radius  as  shown  in 
Equation  (2.53). 
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Figure  2-8:  Geometrical  Relationship  of  Parameters  Used  in  the  Derivation  of  the  Correction 
Factor  for  Open  Sections  in  a  Dome  Structure  (top  view). 
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2  arcsin 


/ 


2  arcsin 


(  ro+Dcos6— a(<$) 
l  cos  0 


)(D- 


( ro+Dcos9— ^ 


cos  9 


a(6) 


(2.53) 


Next,  a (8)  can  be  related  directly  to  the  radius  of  the  hemisphere  (R)  and  the  crush  distance 
of  the  shell  (8)  by  the  following: 


a(6)  =  V#2  -(R-  W  =  \/6(2R  -  S) 


(2.54) 


Finally,  substitution  of  Equation  (2.54)  into  (2.53)  and  the  result  into  (2.51)  yields  the  desired 
result  for  the  correction  factor  (Equation  (2.55)) 


*(«)  = 


1 - arcsin 

7 r 


f  J(^0+D  C°S  e-^/S(2R-6)  ^  (jj  f  r0+Dcos$-y6(2R-6)  \\\ 


cos  9 


y/8(2R  -  6) 


V 


/ 


C  hemisphere 

(2.55) 


Cylindrical  Shell 

The  same  principles  described  in  Section  (2.2.4)  will  also  be  used  to  estimate  the  effect  of 
cutout  sections  in  the  cylindrical  section  with  the  exception  that  the  cutout  shape  will  be  much 
simpler.  The  cutouts  in  the  cylindrical  section  will  be  rectangular  in  shape,  and  thus,  the 
material  removed  will  not  vary  with  axial  crush  distance.  So,  using  Equation  (2.50)  and  again 
assuming  that  £  =  1,  the  correction  factor  of  the  cylinder  can  be  expressed  as 

2-na  —  rupa  „  mp  , 

£  =  - - - =  1  —  —  where  (2.56) 

27 ra  27t 

n  =  Number  of  cutouts 
p  =  Angular  extent  of  cutout 
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Figure  2-9:  Stiffener  Smearing  Process. 

2.2.5  Effect  of  Shell  Stiffeners 

The  energy  dissipation  effects  of  shell  stiffeners  can  be  taken  into  account  using  Equations  (2.31) 
and  (2.49)  by  assuming  that  these  structural  members  are  “smeared”  across  the  surface  of  the 
shell  so  as  to  produce  a  new  effective  shell  thickness.  This  process  is  illustrated  in  Figure  (2-9) . 
If  the  material  properties  of  the  stiffeners  are  identical  to  those  of  the  shell  material,  their  cross 
sectional  area  is  distributed  across  the  surface  of  the  shell  such  that  the  new  effective  thickness 
(teff  )  can  be  determined  by  the  following  expression 

n 

x> 

teff  —  tshell  “I"  27CCL  (2.57) 

where  Ax  is  the  cross  sectional  area  of  the  ith  stiffener,  n  is  the  total  number  of  stiffeners,  a  is 
the  shell  radius,  and  tsheii  is  the  shell  thickness.  If  the  shell  and  stiffener  materials  have  different 
properties,  the  force  distribution  over  the  effective  cross  section  must  be  made  consistent  with 
the  force  distribution  prior  to  the  “smearing”  process  by  using  the  following  relationship 

P ef f  =  P shell  P stiffeners  (2.58) 
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By  substituting  the  appropriate  relationships  and  assuming  that  Peff  will  be  calculated  using 
the  shell  properties,  Equation  (2.58)  can  be  written  as 

n 

27T ateffEshelie  =  2n  at  shellE shelie  +  ^2  A-iEiS 

1=1 

which  can  be  solved  to  determine  the  effective  shell  thickness 

1  71  A 

tejf  =  t shell  +  2-jrnE  "^AiEi  (2.59) 

1  i=l 

When  E\  =  E2  =  ...  =  En  in  Equation  (2.59),  Equation  (2.57)  is  recovered. 
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Chapter  3 


Modeling  and  Analysis  of  ASDS 
Composite  Shell  Structure 


With  the  analysis  tools  developed,  focus  now  shifts  to  the  modeling  and  analysis  phase  of 
the  project.  First,  a  model  will  be  developed  for  the  actual  ASDS  nose  fairing  structure  and 
forward  hemi-head  of  the  pressure  hull  using  simple,  rotationally  symmetric,  geometric  shapes 
for  the  shell.  Stiffeners  and  cutout  sections  will  be  incorporated  using  the  techniques  discussed 
in  Section  (2).  Next,  the  analysis  tools  developed  in  Section  (2)  will  be  applied  to  determine 
the  crushing  force  to  axial  crush  distance  relationship  and  energy  absorption  capability  of  the 
structure.  Finally,  a  maximum  impact  speed  will  be  calculated  by  equating  the  energy  absorbed 
in  the  collision  to  the  kinetic  energy  possessed  by  the  vehicle  just  prior  to  the  impact. 


3.1  ASDS  Modeling 

3.1.1  Actual  Structure 

The  forward  portion  of  the  ASDS  consists  of  a  hemispherical  pressure  hull  and  a  composite 
fairing  structure.  The  hemispherical  pressure  hull  section  is  constructed  of  HY-80  steel  and 
is  welded  to  the  cylindrical  portion  of  the  hull  to  form  a  water  tight  boundary.  The  ASDS 
nose  fairing  (see  Figure  (3-1))  is  a  composite  structure  attached  to  the  forward  end  of  the 
pressure  hull  by  eight  titanium  clevises.  The  structure  has  an  overall  length  of  2.94  m  and  a 
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Figure  3-1:  The  Actual  ASDS  Nose  Fairing  Structure  Showing  Internal  Components  (Viewed 
from  Aft  Looking  Forward). 
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Property 

A 

B 

C 

D 

E 

Ei  (x  106  psi) 

3 

11.6 

4.7 

.0013 

16.4 

E2  (x  106  psi) 

2.7 

10.6 

4.7 

.0013 

16.4 

E3  (x  106  psi) 

.55 

.37 

.45 

- 

- 

^1,2 

.15 

.124 

.15 

.499 

.31 

Fiber  Content 

.5 

.5 

.5 

- 

- 

Tensile  Strength 

Axial  (X)  (ksi) 

39.1 

75 

56 

- 

125 

Transverse(Y)  (ksi) 

39.1 

75 

56 

125 

Compressive  Strength 

Axial  (X)  (ksi) 

33.8 

68 

14 

!  - 

125 

Transverse  (Y)  (ksi) 

33.8 

68 

14 

- 

125 

Shear  Strength  (S)  (ksi) 

9.41 

7.5 

4.5 

- 

77 

A-E-Glass  (1583)/Epoxy  (7780) 

woven  clot! 

1 

B-Carbon  woven  cloth  (6K  5HS  satin 

weave  IM7) 

C- Kevlar  (K285)/  epoxy  (7714)  woven  cloth 

D-Rubber  core  material 

E-Titanium  devise  material  (Ti  6A1-4V) 

Table  3.1:  Properties  of  Materials  Used  in  ASDS  Nose  Fairng. 


pseudo-rectangular  cross  section  with  maximum  dimensions  of  2.21  m  high  and  2.06  m  wide. 
The  structure  itself  consists  of  an  outer  shell  constructed  from  E-Glass  (1583)  /Epoxy  (7780) 
preimpregnated  (prepreg)  woven  cloth  in  a  single  mold  process.  Internally,  the  shell  is  stiffened 
by  a  network  of  longitudinal  and  transverse  woven  carbon  cloth  (6K  5HS  Satin  Weave  IM7) 
beams  which  are  bolted  to  the  shell  and  interconnected  by  titanium  clevis  fittings.  Attached  to 
the  front  end  of  the  shell  is  a  sonar  dome  comprised  of  a  BFGoodrich  Rho-Cor  material  system. 
This  system  consists  of  Kevlar  (K285)/Epoxy  (7714)  prepreg  woven  cloth  with  a  rubber  core 
material.  The  properties  for  the  various  material  systems  described  are  contained  in  Table 
(3.1).  Cutouts  have  been  made  in  shell  and  the  sonar  dome  to  accommodate  such  equipment 
as  thrusters,  the  forward  anchor,  lights  and  video  cameras. 

3.1.2  Simplified  Model 

The  actual  structural  design  of  the  ASDS  had  to  be  simplified  to  rotationally  symmetric 
shapes  before  the  analysis  tools  developed  in  Section  (2)  could  be  applied.  A  representation 
of  the  model  developed  is  shown  in  Figure  (3-2).  For  the  forward  pressure  hull,  no  modeling 
idealizations  needed  to  be  made  since  the  structure  was  already  rotationally  symmetric  and 
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Thruster  Door  Cutotrt- 


.99  R 


l  : 
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-95- 
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— Camera  Cutout 
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.98  R 


Anchor- 

-Transverse  Beams 


-Light  Cutou 


-Forward  Pressure  Hull  Hemi-head 


Figure  3-2:  Model  of  ASDS  Forward  Section  (Dimensions  in  meters). 

hemispherical  in  shape.  The  pressure  hull  has  a  .99  m  radius  and  an  undisclosed  thickness.  The 
actual  thickness  of  the  structure  is  held  as  confidential  by  the  U.S.  Navy  so  as  not  to  reveal 
the  operational  diving  depth  of  the  vessel.  For  this  analysis,  a  notional  thickness  of  19  mm 
will  be  used.  The  sonar  dome  was  modeled  as  a  spherical  dome  section  using  the  axial  depth 
of  the  actual  sonar  dome  and  a  base  radius  equal  to  the  average  of  the  height  and  width  of 
the  actual  sonar  dome.  The  thickness  and  lay-up  construction  was  taken  to  be  that  of  the 
actual  structure.  Three  cutout  sections  simulating  the  locations  of  the  camera  and  lights  were 
located  on  the  model  120°  apart  at  the  same  radial  distance.  The  actual  and  model  sonar 
dome  dimensions  are  summarized  in  Table  (3.2).  The  shell  structure  was  modeled  using  a 
combination  of  a  frustum  for  the  forward  portion  of  the  shell  that  transitions  into  the  sonar 
dome  and  a  stiffened  cylinder  for  the  after  section.  The  radius  of  the  cylindrical  section  was 
taken  to  be  the  average  of  the  maximum  height  and  width  of  the  actual  structure.  The  frustum 
was  used  as  a  transition  section  between  the  sonar  dome  and  the  cylindrical  section  and  uses 
the  corresponding  radius  at  each  end  to  fix  its  dimensions.  For  the  frustum  section,  the  shell 
thickness  was  taken  to  be  that  of  the  actual  material.  The  thickness  of  the  cylindrical  section 
was  taken  to  be  an  effective  thickness  calculated  using  Equation  (2.59)  to  account  for  the  added 
structural  stiffness  of  the  longitudinal  carbon  beams.  The  transverse  beams  were  not  utilized 
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Actual 

Parameter 

Model 

.64 

Axial  Depth  (m) 

.635 

1.91 

Height  (m) 

- 

1.75 

Width  (m) 

- 

1.83 

Average  (m) 

- 

- 

Radius  (m) 

.91 

9.6 

Thickness  (mm) 

9.6 

Camera  Cutout 

.28 

Diameter  (m) 

.28 

Orientation 

0° 

.71 

Radial  Distance  (m) 

.71 

Light  Cutout  1 

.28 

Diameter  (m) 

.28 

-147° 

Orientation 

-120° 

.72 

Radial  Distance 

.71 

Light  Cutout  2 

.28 

Diameter  (m) 

.28 

148° 

Orientation 

120° 

.69 

Radial  Distance 

.71 

Table  3.2:  Comparison  of  Actual  and  Modeled  Sonar  Dome  Dimensions. 

in  the  effective  thickness  calculation  because  their  energy  absorption  in  an  orthogonal  impact 
is  would  be  minimal.  Instead,  these  beams  were  used  to  define  nodal  points  for  buckling  waves. 
Table  (3.3)  summarizes  the  actual  and  modeled  dimensions  of  the  shell  section.  Two  cutout 
sections  1.13  m  high  by  .81  m  wide  were  included  in  the  cylindrical  section  to  account  for  the 
openings  created  by  the  forward  thruster  doors.  (It  is  assumed  that  during  a  collision  imminent 
situation,  the  thruster  doors  would  be  open  and  the  thrusters  in  use.)  Additionally,  the  titanium 
clevises  (used  to  attach  the  internal  beams  together  and  to  the  pressure  hull)  and  the  anchor 
were  treated  as  rigid  objects. 


3.2  Analysis 

3.2.1  Impact  sequence 

As  stated  previously  in  Section  (1),  this  analysis  will  investigate  an  head-on  impact  with  a 
flat,  rigid  wall.  Structural  deformation  will  initiate  at  the  bow  as  the  sonar  dome  is  impacted 
first.  The  deformation  will  progress  .635  m  through  the  sonar  dome,  .95  m  through  the  frustum 
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Actual 

Parameter 

Model 

Cylindrical  Section 

1.36 

Length  (m) 

1.36 

2.21 

Height  (m) 

2.06 

Width  (m) 

2.14 

Average 

- 

Radius  (m) 

1.07 

12.8 

Thickness  (mm) 

35 

Frustral  Section 

.95 

Length  (m) 

.95 

“ 

Fwd  Radius  (m) 
(Sonar  Dome  Radius) 

.91 

- 

Aft  Radius  (m) 
(Cylinder  Section  Radius) 

1.07 

10.1 

Thickness  (mm) 

10.1 

Table  3.3:  Comparison  of  Actual  and  Modeled  Nose  Shell  Section  Parameters. 

section  of  the  nose  shell,  and  .18  m  into  the  cylindrical  section  before  the  anchor  makes  contact 
with  the  wall.  At  this  point,  some  of  the  crushing  load  is  transferred  to  the  pressure  hull  via 
the  anchor  and  simultaneous  deformation  of  the  cylindrical  section  and  the  pressure  hull  results 
for  another  .66  m.  The  deformation  process  ends  when  the  after  most  transverse  beams  and 
titanium  cleaves  make  contact  with  the  wall. 

3.2.2  Sonar  Dome 

Equation  (2.31)  will  be  used  to  determine  the  crushing  force  required  to  deform  the  structure 
throughout  the  progressive  crushing  of  the  sonar  dome.  Implicit  in  the  use  of  this  equation  is 
that  the  structure  can  be  modeled  using  rigid-plastic  techniques.  While  it  is  understood  that  in 
general  composite  structures  behave  in  a  more  brittle  than  ductile  way,  kevlar  composites  have 
been  shown  to  buckle  under  axial  loading  in  modes  similar  to  those  seen  in  metal  structures 
[19].  This  behavior  coupled  with  the  fact  that  the  structure  has  a  rubber  core  was  used  to 
justify  the  use  of  this  analysis  technique  on  this  component. 

The  sonar  dome  is  composed  of  22  plys  of  kevlar  cloth  laid-up  symmetrically  sandwiching  a 
3.81  mm  rubber  core  ([90, 0] 5 , 0 ,  rubber  core,  0,  [90, 0]s)  for  a  total  laminate  thickness  of  9.4  mm. 
The  compressive  and  tensile  flow  stresses  were  calculated  using  a  combination  of  laminated 
plate  theory  and  a  progressive  failure  analysis  [21].  The  detailed  calculations  are  shown  in 
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Appendix  (A)  and,  for  brevity  sake,  only  an  summary  of  the  process  and  the  result  will  be 
included  here.  First,  the  material  properties  of  the  ply  (listed  in  Table  (3.1))  were  used  to 
determine  the  undamaged  0°  ply  stiffness  characteristics.  These  properties  were  then  modified 
using  techniques  of  [21]  to  determine  the  reduced  stiffness  of  the  0°  ply  with  matrix  and  fiber 
damage.  The  stiffnesses  were  calculated  using  the  following  relationship: 


Q  = 


E\  1/7  E\  n 

1— V\V2  1— l/i  1/2  u 

1/ 1  Eo  E7  r\ 

1— V\l/2  \—i/\i/2  U 

0  0  £3 


(3.1) 


The  stiffnesses  were  then  rotated  to  a  90°  ply  orientation.  The  stiffness  of  the  rubber  core 
material  was  also  determined  from  Equation  (3.1)  by  setting  E\  =  E2,  Vi  =  v 2,  and  £3  = 

The  overall  laminate  stiffness  was  then  calculated  using  the  following: 


A* 

Nx 

h 

Q 


1 


(NoohooQoo  +  Ngo°  h$Qc  Q9O0  ~  Nrut,berhTv.bberQrubber)  where 


^laminate 

Total  number  of  “x”  plys 
Thickness  of  individual  ply 
Stiffness  of  individual  ply 


Next,  stress  space  strength  parameters  were  calculated  from  the  material  test  data  provided  in 
Table  (3.1)  using  Equations  (3.2)  and  (3.3)  and  the  generalized  von  Mises  value  of  —5  for  jF\*2. 
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Fi 
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■PT.2 
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(3.2) 


(3.3) 


The  strain  space  strength  parameters  and  failure  envelopes  could  then  be  determined  through 
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the  relationships  in  Equations  (3.4)-(3.6). 


Gk,i 

—  FijQifiQjj 

(3.4) 

Gj 

=  FiQij 

(3.5) 

Gitj£iSj  1  Gj£i 

=  1 

(3.6) 

A  progressive  failure  model  [22]  was  then  used  to  approximate  the  stress  strain  relationship  of 
the  material  under  pure  compressive  and  tensile  loadings  (oriented  in  the  0°  ply  direction).  The 
progressive  failure  model  allows  for  a  rough  prediction  of  post  first  ply  failure  behavior  of  the 
laminate.  The  model  works  by  calculating  the  strain  of  first  ply  failure  for  the  assumed  load 
condition  and  then  determines  the  stress  using  the  undamaged  laminate  stiffness  ({<r}  =  A*{c}). 
After  first  ply  failure,  the  limiting  ply  is  failed  to  a  degraded  state,  a  new  laminate  stiffness  is 
calculated,  and  the  next  ply  failure  strain  and  stress  is  determined.  A  ply  can  be  damaged  at 
most  twice  -  once  by  degraded  matrix  and  the  final  by  fiber  failure.  In  certain  cases,  based  on 
the  properties  of  the  plys  used,  the  degraded  ply  failure  strains  can  be  less  than  that  of  the 
undamaged  ply.  For  this  condition,  the  ply  failure  occurs  after  its  undamaged  failure  envelope 
is  surpassed.  The  tensile  and  compressive  stress  versus,  strain  estimates  are  shown  in  Figure 
(3-3). 

The  compressive  (<jco)  and  tensile  ( ato )  flow  stresses  were  then  determined  to  be  8.41  MPa 
and  53.35  MPa  respectively  using  Equation  (3.7). 


Now,  the  crushing  force  (PSOnar  dome)  and  cumulative  energy  absorption  (EsonaT  dome)  of  the 
sonar  dome  can  be  determined  from  Equation  (2.31)  and  (3.8)  with  the  cutout  correction  factor 
(Equation  (2.55))  applied  to  the  crushing  force  over  the  appropriate  regions  of  the  shell. 

E=  f6  P(6)d6  (3.8) 

Jo 

The  resulting  crushing  force  and  cumulative  energy  absorption  profiles  for  an  axial  crush  dis¬ 
tance  (6)  are  shown  in  Figure  (3-4).  Here  it  is  seen  that  the  total  energy  absorbed  for  a  crush 


1 

<7o  =  — 

Ef  JO 
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300.000 


Figure  3-3:  Compressive  and  Tensile  Stress  vs  Strain  Estimates  Calculated  using  the  Progressive 
Failure  Method. 

distance  of  .635  m  is  .058  MJ.  This  figure  also  clearly  shows  the  effect  of  the  cutouts  on  the 
crushing  force  curve. 

3.2.3  Frustum  Section  of  Nose  Shell 

The  frustum  section  is  composed  of  a  22  ply,  E-Glass  (1583)/Epoxy  (7780),  woven  cloth  com¬ 
posite  system  laid  up  symmetrically  ([90,45,0,  — 45)2,90,45,0,0,45,90,  [90,45,0, —45)2)  for  a 
laminate  thickness  of  10.1  mm.  The  analysis  of  the  frustum  section  was  done  by  approximating 
the  frustum  as  a  cylinder  with  a  radius  equal  to  the  average  of  the  forward  and  after  radii  of  the 
frustum  and  using  the  analysis  discussed  in  Section  (2.2.3).  The  details  of  these  calculations  are 
contained  in  Appendix  (B).  In  order  to  use  this  technique,  the  appropriate  circumferential  (Eg) 
and  flexural  (E/)  Young’s  moduli  needed  to  be  determined.  The  method  described  in  Section 
(3.2.2)  was  initially  used  to  determine  the  strain  space  failure  envelopes  for  each  of  the  ply  ori¬ 
entations  used  in  the  laminate.  Subsequently,  the  failure  envelopes  were  used  to  develop  an  £2 
versus,  tensile,  transverse  strain  relationship  using  the  progressive  failure  technique  described 
previously.  £2  was  chosen  as  the  circumferential  strain  Young’s  moduli  because  as  the  buckling 
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0.08  - 1 - i -  0.16 

Effect  of  Cutout  Sections 


Figure  3-4:  Crushing  Force  and  Cumulative  Energy  Absorption  as  a  Function  of  Axial  Dis¬ 
placement  for  the  Sonar  Dome. 

lobes  are  formed  and  bulge  outward  and  inward,  the  laminate  will  be  strained  in  the  2-direction 
with  respect  to  the  laminate  orientation.  In  reality,  both  compressive  (bulging  inward)  and 
tensile  (bulging  outward)  strains  will  be  experienced  and  thus,  different  E<i  versus,  strain  re¬ 
lationships  will  result.  However,  because  the  compressive  and  tensile  strengths  of  the  plys  are 
approximately  equal  (Table  (3.1)),  the  difference  is  small  and  only  the  tensile  Eo  versus,  strain 
relationships  will  be  used  to  simplify  the  analysis.  Figure  (3-5)  shows  the  Eo  versus,  strain 
relationship  that  was  determined  as  well  as  the  order  of  ply  failure. 

The  actual  circumferential  strain  was  then  calculated  as  a  function  of  x  and  6  using  Equa¬ 
tions  (2.36),  (2.48),  and  (2.44).  This  relationship  is  shown  in  Equation  (3.9). 


x 

a H ' 


x 6  jy/at 


ay/at 


1 

4 


(3.9) 


Since  the  circumferential  strain  varies  both  with  x  and  6  and  £2  varies  with  S,  Eg  will  vary 
with  both  x  and  6.  Eg  was  then  calculated  for  discrete  values  of  x  and  6  in  preparation  for  a 
numerical  integration  scheme  to  be  used  in  Equation  (2.42). 

To  determine  E /,  a  relationship  between  flexural  strain  (s/)  and  the  axial  crushing  distance 
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Transverse  Strain  (in/in) 


Figure  3-5:  Transverse  Young’s  Modulus  Variations  with  Pure  Tensile.  Transverse  Strain  for 
Frustum  Section. 

(<5)  was  first  developed  using  Equations  (2.38),  (2.39),  and  (2.48).  (see  Equation  (3.10)) 

(3.10) 

Equation  (3.10)  was  then  used  to  determine  a  ej  versus.  6  curve  for  each  ply  over  the  entire 
crushing  distance  for  one  buckling  lobe  (2 Vat).  For  each,  the  strain  was  determined  at  the 
midplane  of  the  individual  plys.  Using  these  plots  and  the  strain  failure  envelopes,  the  dis¬ 
placement  ait  which  each  ply  failed  in  tension  could  be  determined.  Figure  (3-6)  shows  how  this 
process  was  done  for  the  0°  oriented  plys  in  the  laminate. 

Flexural  stiffness  (D*),  E /,  and  laminate  thickness  (t)  were  then  recalculated  after  each  ply 
failure  using  Equations  (3.11)-(3.13). 


m  =  p  f)iQ}z2dz 

(3.11) 

E,  =  {(l^T1),,,}'’ 

(3.12) 

t  =  nh  where 

(3.13) 
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Figure  3-6:  Flexural  Strain  to  Axial  Crush  Distance  Relationships  for  the  0  Degree  Oriented 
Plys  in  the  Frustum  Laminate. 

n  =  number  of  intact  plys 
h  =  ply  thickness 

Figure  (3-7)  shows  the  variations  in  E /  and  t  as  6  increases.  Notice  also  that  as  more  plys  fail, 
Ef  becomes  more  erratic  and  diverges  further  and  further  from  the  in-plane  modulus.  This 
is  expected  since  the  laminate  starts  out  homogenized  and  gradually  de-homogenizes  as  each 
successive  ply  fails. 

The  crushing  load  and  cumulative  energy  absorption  can  now  be  determined  as  a  function 
of  axial  crushing  distance  (Figure  (3-8))  using  Equations  (2.42)  and  (3.8).  Since  neither  E$ 
nor  Ef  was  constant  or  equal  during  the  deformation  process,  the  simplified  form  of  Equation 
(2.42)  could  not  be  used  here.  Figure  (3-8)  shows  that  five  buckling  lobes  form  as  the  frustum 
is  crushed  and  that  a  total  energy  absorption  of  .105  MJ  is  realized. 
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Figure  3-7:  Flexural  Young’s  Modulus  and  Laminate  Thickness  Variation  with  Axial  Crush 
Distance  of  the  Frustum  Section. 
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Figure  3-8:  Crushing  Load  and  Cumulative  Energy  Absorption  Characteristics  of  the  Frustum 
Section. 
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3.2.4  Cylindrical  Section  of  Nose  Shell 

This  outer  shell  of  the  stiffened,  cylindrical  section  is  composed  a  28  ply,  E-Glass  (1583) /Epoxy 
(7780),  woven  cloth  composite  system  laid  up  symmetrically 
([90,45,0, -45)2,90,45,0,0,45,90,  [90,45,0, —45)2)  for  a  laminate  thickness  of  12.8  mm.  Ad¬ 
ditional  support  is  provided  by  an  internal,  frame  structure  constructed  using  carbon,  woven 
cloth  (6K  5HS  satin  weave  IM7),  composite,  C-channel  beams.  The  beams  are  32  plys  thick 
for  a  laminate  thickness  of  11  mm  and  laid-up  symmetrically  in  the  following  configuration: 
[90, 45, 0,-45]4[90, 45, 0,-45]4.  The  shell  analysis  was  accomplished  by  the  same  process  as  de¬ 
scribed  in  Section  (3.2.3)  and  is  shown  in  detail  in  Appendix  (C).  However,  the  effect  of  the 
stiffeners  still  needed  to  be  taken  into  account.  This  was  done  by  "smearing”  the  stiffeners 
into  the  shell  and  then  calculating  an  effective  thickness  (te//)  using  Equation  (2.57).  The 
detailed  stiffener  and  t eff  calculations  are  included  in  Appendices  (C)  and  (D)  with  the  t eff 
determined  to  be  35  mm.  An  effective  buckling  lobe  length  was  then  determined  based  on  t eff. 
The  crushing  force  and  cumulative  energy  absorption  could  then  be  determined  using  Equations 
(2.42)  and  (3.8).  Figure  (3-9)  shows  the  force  and  energy  profiles  at  various  displacements  and 
also,  that  the  cylinder  was  only  allowed  to  crush  .84  m  in  two  buckling  lobes  while  Table  (3.3) 
indicates  that  the  section  is  1.36  m  long  and  should  buckle  in  three  lobes.  This  reduction  in 
length  was  due  to  the  effects  of  the  transverse  beam  structures.  The  forward  beam  and  titanium 
fittings  served  to  force  a  nodal  point  at  the  beginning  of  the  thruster  door  section  as  seen  in 
Figure  (3-2).  The  after  beams  and  rigid  titanium  clevises  at  the  aft  portion  of  the  thruster  door 
create  another  node  point  and  serve  to  terminate  the  deformation  process.  So,  only  the  .84  m 
shell  and  stiffener  portion  between  the  forward  and  after  beams  will  deform  during  the  impact. 
.49  MJ  of  energy  are  absorbed  during  this  deformation. 

3.2.5  Analysis  of  Forward  Hemi-head  of  Pressure  Hull 

The  forward  pressure  hull  hemi-head  is  a  metallic  structure  constructed  from  HY-80.  As  pre¬ 
viously  stated,  the  thickness  used  in  the  analysis  will  be  assumed  to  be  19  mm  due  to  the 
confidentiality  of  the  actual  thickness.  The  deformation  starts  as  the  forward  anchor  contacts 
the  collision  boundary  and  is  pushed  back  into  the  pressure  hull.  The  deformation  is  assumed 
to  be  of  the  form  described  in  Section  (2.1)  and,  as  such,  Equation  (2.31)  will  be  used  to  de- 
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Figure  3-9:  Crushing  Force  and  Cumulative  Energy  Absorption  of  the  Cylinder  Section  of  the 
Nose  Shell. 

termine  the  crushing  force.  The  flow  stress  was  calculated  from  a  partial  stress  versus,  strain 
curve  using  Equation  (3.7).  The  detailed  calculations  and  the  stress  versus,  strain  profile  used 
are  included  in  Appendix  (E),  and  the  crushing  force  and  cumulative  energy  absorption  versus, 
displacement  results  are  shown  in  Figure  (3-10).  Throughout  this  analysis,  it  is  assumed  that 
the  forward  hemi-head  remains  water  tight.  Further  analysis  must  be  performed  to  determine 
the  point  at  which  the  pressure  hull  fails.  By  allowing  the  pressure  hull  to  deform  to  the  full 
axial  extent  allowed  for  the  nose  shell,  7.566  MJ  of  energy  will  be  absorbed. 

3.2.6  Kinetic  Energy  Calculations 

After  determining  the  energy  absorption  capability  of  each  modeled  nose  fairing  component,  a 
determination  of  the  kinetic  energy  possessed  by  the  ASDS  was  needed  for  various  speeds.  The 
kinetic  energy  was  determined  using  Equation  (3.14) 

KEaSDS  ~~~  2  (mvehicle  "b  added) V  (3.14) 

where  is  the  submerged  mass,  matided  is  the  added  mass,  and  v  is  the  speed  of  the 

vehicle. 
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Figure  3-10:  Crushing  Force  and  Cumulative  Energy  Absorption  for  the  Forward  Pressure  Hull 
Hemi-head. 
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Figure  3-11:  Description  of  Standard  Submarine  Weight  Terms 


Parameter 

Weight 
(metric  tons) 

NSC 

56.2 

Main  Ballast  Tanks 

4.8 

Free  Flood 

27 

Aenv 

88 

Table  3.4:  ASDS  Weights 


The  submerged  mass  of  the  vehicle  was  calculated  using  the  ASDS  weight  report  [16]. 
Weight  reports  for  submarines  are  standardized  and  broken  down  into  categories  of  weights 
based  on  the  operating  condition  of  the  submarine.  Figure  (3-11)  shows  a  brief  summary  of 
the  standard  terminology  used  in  submarine  weight  analysis  and  how  each  relates  to  the  other. 
Condition  A  is  the  dry  vehicle  weight.  When  the  variable  load  (variable  ballast  tank  water, 
personnel,  etc.)  is  added,  the  new  weight  created  is  the  normal  surfaced  condition  (NSC).  The 
main  ballast  tanks  are  filled  with  water  to  submerge  the  vehicle;  so,  the  weight  of  this  water 
added  to  the  NSC  becomes  the  submerged  displacement  (Asub).  Finally,  during  submergence, 
non-water  tight  free  flood  areas  (i.e.-between  the  fairings  and  the  pressure  hull)  fill  with  water. 
This  water,  while  in  pressure  equilibrium  with  the  surrounding  water,  moves  with  the  submarine 
as  it  transits,  and  when  its  weight  is  added  to  A sub  the  resultant  is  the  envelope  displacement 
( Aenv ).  For  this  case,  it  is  this  Aenv  that  is  required  for  the  analysis.  Table  (3.4)  shows  the 
various  weights  used  to  determine  the  envelope  displacement. 

Next,  the  added  mass  was  determined.  Since  this  vehicle  is  relatively  slender,  it  is  expected 
that  the  added  mass  in  surge  will  not  be  large.  However,  two  approximations  were  made  using 
a  sphere  and  an  ellipsoid.  The  size  of  the  sphere  was  set  such  that  the  frontal  area  was  equal 
to  that  of  the  ASDS  (Figure  (3-12)).  Once  the  sphere  dimensions  are  known,  the  added  mass 
can  be  calculated  using  Equation  (3.15) 


^addcd  sphere  — 
(where  pwater  = 


t^Pwater^sphere  Pwater  ^-sphere  ^.5  metric  tons 

1.025 metric,t0nS) 
m6 


(3.15) 


The  added  mass  of  the  vehicle  was  also  approximated  using  an  ellipsoid  with  the  frontal  area 
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Figure  3-12:  Dimensions  Used  During  Added  Mass  Calculations 

equal  to  that  of  the  ASDS  (as  in  the  case  of  the  sphere)  and  a  major  axis  length  equal  to  that 
of  the  overall  vehicle  (Figure  (3-12)).  The  added  mass  of  the  ellipsoid  can  then  be  determined 
using  Equation  (3.16)  and  the  graphs  from  [17]  for  the  added  mass  coefficient  (Ctt) 

4 

Padded  ellipsoid  =  Ca-irpwaterab2  =  2.1  metric  tons  (where  Ca  =  .13)  (3.16) 

Since  the  actual  vehicle  shape  is  not  as  bluff  as  a  sphere,  yet  not  as  fine  as  an  ellipse,  the  added 
mass  of  the  vehicle  lies  somewhere  between  the  two  calculated  values.  So.  the  average  value 
was  used  as  the  vehicle  added  mass.  Now,  the  kinetic  energy  possessed  by  the  ASDS  at  various 
speeds  can  be  determined  using  Equation  (3.14)  as  shown  in  Figure  (3-13). 

3.3  Results 

With  the  energy  absorption  of  each  section  determined,  the  overall  energy  absorption  of  the  nose 
structure  can  be  determined  by  summing  the  energy  absorbed  by  each  section  throughout  the 
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Figure  3-13:  Vehicle  Kinetic  Energy  for  Various  Speeds. 

impact  sequence.  Figures  (3-14)  and  (3-15)  show  the  integrated  crushing  force  and  cumulative 
energy  absorption  profiles  for  the  entire  deformation  process.  Two  curves  are  shown  in  Figure 
(3-15)  depicting  the  energy  absorption  of  the  nose  structure  with  and  without  the  effects  of  the 
pressure  hull  included.  Additionally,  the  kinetic  energy  curve  shown  in  Figure  (3-15)  is  overlaid 
on  this  plot.  Using  the  curves  in  Figure  (3-15)  initial  speed  versus,  crush  distance  assessments 
can  be  made.  Without  the  pressure  hull  included,  the  nose  structure  can  absorb  .67  MJ  of 
energy  which  would  allow  the  complete  absorption  of  all  kinetic  energy  possessed  by  the  vehicle 
with  an  initial  speed  of  7.5  knots.  If  the  effects  of  the  pressure  hull  are  included,  a  total  of 
8.235  MJ  of  energy  can  be  absorbed.  This  equate  to  an  initial  vehicle  velocity  of  26.4  knots. 
As  previously  stated,  the  point  at  which  the  pressure  hull  will  fail  during  its  deformation  was 
not  determined;  so,  a  the  safe  impact  speed  for  the  vehicle  cannot  yet  be  determined.  However, 
when  this  point  is  determined,  Figure  (3-15)  will  provide  this  safe  speed  estimate. 
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Figure  3-14:  Crushing  Force  Profile  for  the  Entire  Deformation  Process. 
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Figure  3-15:  Total  Cumulative  Energy  Absorption  and  Vehicle  Kinetic  Energy. 
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Chapter  4 


Summary  and  Conclusions 


Throughout  the  course  of  this  analysis,  several  new  deformation  models  have  been  developed. 
The  force  required  to  deform  a  hemispherical  shell  loaded  axially  with  a  rigid  plate  by  axisym- 
metric  dishing  was  related  to  shell  thickness,  material  flow  stress  (compressive  and  tensile), 
and  axial  deflection  using  a  rigid-plastic  analysis  technique.  Additionally,  a  deformation  model 
for  a  cylindrical  shell  composed  of  elastic  material,  crushed  axially  by  a  rigid  flat  plate,  was 
proposed.  This  model  assumed  a  concertina  type  deformation  mode  and  used  shell  thickness, 
cylinder  radius,  material  stiffness,  and  axial  displacement  to  determine  the  force  required  to 
crush  the  shell.  Further,  methods  were  set  forth  to  include  the  effects  of  cutout  sections  and 
internal  stiffening  members  on  the  crushing  force  characteristics  of  the  shells.  The  effect  of  the 
cutout  section  was  taken  to  be  a  reduction  of  the  intact  shell  crushing  force  by  a  ratio  equal  to 
the  shell  area  remaining  divided  by  the  total  area  of  the  intact  shell  in  the  local  vicinity  of  the 
cutout.  The  effect  of  internal  stiffeners  was  incorporated  by  increasing  the  thickness  of  the  shell 
material  to  a  new  effective  thickness  determined  by  making  the  load  carried  by  the  effective 
structure  equal  to  the  combined  load  on  the  stiffeners  and  shell. 

These  tools  were  used  to  determine  the  crushing  force  and  energy  absorption  characteristics 
of  the  composite  nose  fairing  structure  and  forward  pressure  hull  of  the  U.S.  Navy’s  Advanced 
Seal  Delivery  System  in  a  head-on  impact  with  a  rigid  wall.  The  complex  structure  of  the 
ASDS  was  first  simplified  and  modeled  to  permit  the  use  of  the  tools  described  above.  The 
shell  structures  were  approximated  as  rotationally  symmetric  shells  (i.e.-  cylinders,  frustum, 
hemispherical  domes)  with  the  internal  shell  stiffeners  being  incorporated  into  the  model  in 
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two  ways.  The  energy  absorption  of  the  longitudinal'elements  was  captured  by  increasing  the 
effective  thickness  of  the  shell  in  the  vicinity  of  the  longitudinal  beams,  while  the  transverse 
stiffening  elements  served  only  to  force  nodal  points  in  the  shell  buckling  pattern  since  their 
energy  absorption  during  a  head-on  collision  is  minimal.  Further,  major  internal  components, 
such  as  the  titanium  clevises  (used  to  connect  the  internal  beams)  and  the  forward  anchor,  were 
treated  as  rigid,  massless  objects. 

The  structure  was  divided  into  a  hemispherical  sonar  dome,  a  frustum  and  stiffened  cylin¬ 
drical  nose  shell,  and  a  hemispherical  forward  pressure  hull  based  on  the  shape,  material  com¬ 
position,  and  thickness  of  the  actual  structure.  Using  laminated  plate  theory  and  a  progressive 
failure  model,  stiffness  versus  crush  distance  and  material  thickness  versus  crush  distance  rela¬ 
tionships  were  developed  for  each  of  the  material  systems  used.  These  relationships  were  then 
used  in  conjunction  with  the  rigid-plastic  hemispherical  and  elastic  cylindrical  shell  models 
to  predict  the  force  deflection  and  energy  absorption  of  each  component.  The  component- 
level  energy  absorptions  were  then  summed  to  approximate  these  characteristics  for  the  overall 
structure.  Finally,  the  absorbed  energy  was  equated  to  the  vehicle's  pre-impact  kinetic  energy 
in  order  to  determine  an  initial  speed  just  prior  to  the  collision.  The  composite  nose  fairing 
structure  was  deemed  capable  of  absorbing  the  impact  energy  associated  with  a  vehicle  speed 
of  7.5  knots.  If  the  effects  of  the  pressure  hull  were  included,  the  structure  could  absorb  the 
kinetic  energy  from  a  26.4  knot  impact. 

4.1  Areas  Additional  for  Study 

This  analysis  has  taken  the  first  step  in  determining  the  energy  absorption  capability  of  the 
ASDS.  However,  considerable  work  remains  to  be  done  before  a  final  assessment  can  be  made 
of  the  vehicle’s  overall  crashworthiness.  The  following  is  a  list  of  topics  for  further  study  in  this 
project: 

•  Model  validation  and  calibration:  The  models  developed  herein  have  not  been  cali¬ 
brated  or  confirmed  via  actual  test  data.  Cylindrical  and  hemispherical  composite  sections 
of  the  appropriate  length  to  radius  and  thickness  to  radius  ratios  need  to  be  constructed 
for  each  composite  material  and  crushed  under  quasi-static  and  dynamic  conditions.  The 


55 


results  should  be  used  to  confirm,  improve,  and  calibrate  the  methods  developed  in  this 
thesis. 

•  Rupture  Point  of  Forward  Pressure  Hull:  The  overall  objective  of  the  project  is  to 
determine  a  safe  impact  speed  for  the  vehicle  such  that  during  an  impact  the  pressure  hull 
will  not  rupture  and  the  crew  will  survive  the  collision.  The  rupture  or  collapse  point  of 
the  forward  pressure  hull  needs  to  be  determined  and  may  vary  with  depth  of  the  vehicle. 

•  Strength  Analysis  of  Nose  Fairing  Mounting  Plates:  The  composite  nose  fairing 
is  mounted  to  reinforced  metal  plates  that  are  welded  to  the  forward  hemi-head  of  the 
pressure  hull.  An  analysis  needs  to  be  done  to  determine  whether  these  structures  will  be  / 
capable  of  withstanding  the  force  of  the  impact  or  if  they  will  bend  or  shear  off  completely. 

•  Energy  Absorption  of  Internal  Components:  The  energy  absorption  potential  of  the 
internal  components  (such  as  clevises,  angle  bulkheads,  the  anchor,  the  sonar  equipment, 
etc.)  have  not  been  included  in  this  analysis.  The  effects  of  each  of  these  components 
can  be  determined  individually  and  added  to  the  curves  developed  in  this  report  to  get  a 
more  accurate  approximation  of  the  energy  absorption  of  the  vehicle. 

•  Effect  of  Impact  Angle  on  Energy  Absorption:  This  analysis  looked  at  an  “impact 
sequence”  for  a  head-on  impact  angle.  In  reality,  this  situation  is  likely  never  to  happen. 

A  relationship  should  be  developed  for  various  impact  angles  using  model  testing  or  other 
means  to  relate  the  reduction  in  energy  absorption  that  is  realized  by  an  off-axis  impact 
to  that  of  an  orthogonal  impact. 
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Appendix  A 


Sonar  Dome  Analysis 
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Ply  Stiffness  and  Layup  Inputs: 


1.  Material  System:  Kevlar  (K285)/  epoxy  (7714)  woven  cloth 

2.  Layup  Configuration:  [90,0]5,0101[90,0]5 


n0;=I2  Nn45:=o 

N  90-10  Np45:=0 

3.  Total  Number  of  Plys  Totai_piys:=  n  0+n  n45+N  ^5+N  90 

4.  Ply  Thickness  h—.Olin  h  rubber:=  15-in 

5.  Laminate  Thickness  t  L  :=Totai_piysh+h  ^ber 

6.  Strain  Applied  to  the  Material  (circumferential): 

P=atan(ey/ex)  p  :=o-deg 

Shear  Strain:  e ,  :=o 


7.  Number  and  Orientation  of  Ply  Angles: 

Plys  -4 
a:=l..Plys 
c  !=1  ..3*Plys 


e  := 


‘  0 

90 

45 

-45 


•deg 


8.  Engineering  Constants: 
Kevlar 


E  x  :=4.7-MSI 


Ey  :=4.7-MSI 


vf:=5 


E  m  :=3.4  GPa 


Gm:=1.26GPa 


Rubber 

E  rubber  ■  ^  304  *psi 

v  rubber”^ 


1  L  =  9.398°mm 
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Model  Parameters 


1 .  Radius  of  Base  of  Sonar  Dome: 

2.  Axial  Depth: 

3.  Radius  of  Curvature 

4.  Radius  to  Center  of  Cutouts 

5.  Axial  Distance  to  Center  of  Cutout 

6.  Slope  of  Dome  Surface  at  Cutout  Center 

7.  Diameter  of  Cutout 

8.  Radius  where  Cutout  Starts 

9.  Radius  where  Cutout  Ends 
10  Number  of  Cutouts 


R  model  ,-36  in 
®  model 25in 

2  2 
._!  D  model  model 

K  curvature*  7' - r - — - 

L  u  model 


R  curvature  “  0.976^m 


R  cutout*  27  in 


e“^R  curvature^  “  R  cutout^  D  cutout  ~  0.282*m 


“D  cutout +R  curvature 


cutout  cutout curvature 


0)  =  58,003°deg 


d  cutout 

r0:=23.693in 
r  f  t=29.66-in 


n  cutouts 

Ply  Strength  Parameter  Inputs: 


1.  Tensile  Strength  in  Longitudinal  Direction 

2.  Compressive  Strength  in  Longitudinal  Direction 

3.  Tensile  Strength  in  Transverse  Direction 

4.  Compressive  Strength  in  the  Transverse  Direction 

5.  Shear  Strength 


x:=56-ksi 


X  prime :='4  KSI 


Y:=56KSI 


Y  prime  :=14KSI 


S  :=4.5-KSI 


Matrix  and  Fiber  Damage  Inputs 

Ef  star  *1  y  **” ’5161  n*—.l 
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Damage  Ply  Calculations: 


E  fx  • 


=  Ex 

"vf 


Gfx:“ 


1  .  Vjn 
l  -nn  s’ — 

5  V  f 

v  ml 
-ns“ 

Es 

^  m 

•H 


Efy‘=l 


i-l 


l+ny  Vp 


E  y_damagerfE) 


1 

v  m  ] 

T1  V' - 

1  1  y  vf 

i  .  V  m 

‘-MyVF 

E  fy  F  E  m  j 

E  s_damaged(E)  •“ 


i+*H  s  — 


_L+JL 

v  m 

T\  s- - 

S  Vf 

Gfx  +  F 

G  m 

rl 


Building  Ply  Data: 

Undamaged  ply,  ply  with  m 


Undamaged  ply,  ply  with  matrix  damage  and  ply  with  fiber  damage. 


E  xx  • 

'  Ex  ' 
Ex 

Eyy 

Ey 

E  y_damaged  E  m_star 

E  ss  := 

E  s 

E  s_damaged  E  m_star 

E  X-E  f_star  _ 

E  y_damaged  E  f_star 

_  E  s_damaged  E  f_star 

V  X 

E  xy_star 

v  xx:= 

v  x'E  m_star 

EE  xy_star  * 

E  xy_star'E  m_star 

_  v  x’E  f_star  _ 

E  xy_star'E  f_star 

xxq:=x 


XX primeq  ,-X prime’  -fe" 


Wq  :=Y 


YY  primeq  ,-Y  prime 


ssq  :=s 


v  y>q  ,-v  xx, 


cyyq 

q  Exxn 


Rotation  Matrices: 


mla:=cos6a  '  na  '=s‘n  ®a 


J  := 


ml  2  n  2 
a  a 

n  2  ml  2 


2-mI  *n 
a  a 


-2mlana 


•ml  n  ml.-n.  ml  2-  n„  2 
a  a  a  a  a  a 


Rotation2a  !=| 


ml 


ml 


ml  -no  -mi  n 
a  a  a  a 


Rotation  1  != 


ml 


ml 


.2,2  .  2  „  2 
ml  •  n  ml  •  n 
a  a  a  a 


mIa2"a2  m,a2na2 


2-  ml  2-  n  2 
a  a 

2  ml  2-  n  2 
a  a 


-2  ml,  2-  na  2 


4  ml  2*  n  2 
a  a 

4  mla  2'  na  2 
-4  ml  2-  n  2 


ml 


ml  •  n 
a  a 


[  mla  2-  na  2] 

na  -mla'  na  3  mla'  na  3~  m,a  3’na  2[mla'  na  3“  mla  3"a] 

[  “I,  3’na-mla-  na  3] 


mla  3’na 


mla  3'na— mIa'  na  3  2f  mla  3  na-mla' 
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Ply  Stiffness  Calculation: 

Calculates  the  zero  direction  stiffness  matrix  for  the  various  ply  states  (undamaged,  matrix  damaged,  and 
fiber  damaged) 


Exxq  vy>’qExxq  ^  1 

^  rubber  v  rubber'^  rubber  n 

“  T  !  F“ 

1  —  v  rubber  1  —  v  rubber 

v  rubber'^  rubber  E  rubber  n 

1  ~v  xxq'v  y>’q  1-v!Oiq'vyyq 

vxXq'^yyq  Eyyq  Q  rubber'- 

2  ^  ^ 

1  —  v  rubber  1  —  v  rubber 

n  n  2  E  rubber 

|-vxxqvyyq  1-vxxqvyyq 

0  0  EssqJ 

1  rubber 

Creates  the  stiffness  matrices  for  every  combination  of  ply  angle  and  ply  state.  Qnew  is  a  3  x  4  matrix 


!■ -Rotation  lfl 


1,1 

2,2 

1,2 

3,3 


Qrotq,a  i 

^rotq,a  3 

^  rotq,a  5 

Qr0‘q,a  3 

Qrotq,a  2 

Q  rotq,a  6 

Q  rotq,a  5 

Qro’q,a6 

^  rotq,a  4 

Stress  Space  Strength  Parameters: 

A  separate  F  must  be  calculated  for  each  ply  state 


FV=' 


FF 


xy_stara 


xxq-xxprimeq 


FF 


xv  starn 


^XXq.XXpnmeq.YYqTYprjmeq 


1 


^XXq.XXprimVYYq.YYprim^ 


YY  -YY 


q Pnmeq 
0 


SS„ 


FV=I 


1  _  1 

XXq  XX  primeq 

1  _  1 
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Rotation  12 

1  2,2 
G  basel i  2 

G  baselj  3  3 
Gbasel2  j  j 

Rotation^ 

G  basel2  2  2 
G basel 2  j  2 

Gbasel2  3  3 

G  basel 3  j  j 

Rotation^ 

G  basel 3  ^  2 

G  basel 3  ]2 

G  ltrans*” 

G  basel 3  3  3 

G  basel]  j  j 

Rotation  I-j- 

G  basel]  2  2 
Gtose1,  ,  2 

G  basel]  3  3 

G  basel 2  j  j 

Rotationly 

G  basel 2  2  2 
Gbasel2  j  2 

Gbasel2  3  3 
Gbasel3  ^  ( 

Rotation  ly 

°basel3  2  2 

G  basel 3  j  2 

G  basel 3  3  3 

G  basel j'j  j 

Rotation  1^* 

G  basel]  2  2 

G  basel]  j  2 

G basel]  3  3 
Gbasel2  j  j 

Rotation^- 

G  basel 2  2  2 
Gbasel2  j  2 

Gbasel2  3  3 

Gbasel3  ,  j' 

Rotation^- 

G  basel3  2  ? 
Gbasel3  {  2 

G  basel3  3  3 

Rotation^ 

i 

G  base23  2 

Rotation24* 

Gbase2]  j 
Gbase2]  2 

Rotation24- 

G  base22  j 

G  base22  9 

Rotation24- 

G  base23  j 

G  base23  2 
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G  ltransj. 

®  ltransj,  ^ 

G  Itrans^  5 

G  2trans^ 

II 

G  ltran%  ^ 

G  ltransc  ^ 

Gltranstg 

^tran^ 

G  ltran^  ^ 

G  ltran^  ^ 

0  ltran^-  4 

®  20^^ 

Formation  of  Strain  Space  Failure  Envelopes: 


1 .  The  following  relationship  is  used  to  determine  the  failure  envelope:  g; 

2.  After  shifting  to  polar  coordinates,  the  above  equation  can  be  solved  for  r  using  the  quadratic  equation: 

£j®r-cos$  E2*r-sin<fr  Ar^-j-B-r— 1=0 

3.  Calculation  of  coefficients  for  the  quadratic  formula 

A(<t>  ,b)  :=  Gb  2-sin(i)>)2+2-  Gb  |  ^ -cos( ^ ) -sin( * )  +-  Gb  ]  t  -cosC ()>  )2 
B  4> , b.e  s  :=2-  Gb  ^  3 -sin( * )  e  s-f-  Gb  ^  4-sin(<|>)+2-  Gfe  (  3<os(*)-s  s+-  Gfe  ^  4«os(4») 

C  4>  ,b,e  s  Gb  ^  4*  s-(-  0|,  3  s  -  I 


r  , b, e  s 


B  $,b,es  +aJb  4»,b,es  4A(<|>,b)C  4>,b,Es 


Strain  Space  Failure  Envelopes: 

<t>  .'=0-deg,l-deg..360deg  y(x)  :=tan(p)-x  x  :=0,co^^  ,.cos(  P)  sx4>,b,£s  :=r  ,b,c  s  -cos(  4> )  £y+,b,ss  :=r^,btEs  sin(  <t> ) 


Ply  Strain  Strength  Envelopes 


A  0 

3— 

( 

t 

y  \  i  i  i  i 

-0.04  -Tt :0j - - Oof — 

0^ - ^jQ&f 

- .  A  A’ 

v.U; 

_ 

:  j 

Longitudinal  Strain 


-  0  ply-Undamaged 

.  0  ply-matrix  damaged 

.  90  plv-Undamaged 

-  90  plv-matrix  damaged 

-  45  ply-Undamaged 
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Calculation  of  Degraded  Laminate  Stiffness: 

The  first  ply  failure  strain  is  first  calculated.  Subsequently,  a  progressive  failure  of  plys  is  allowed  and  the 
stiffness  calculated  based  on  the  ply  damage  history. 


pioti  := 


al*-l 

bl«-I 


count*-  2 
nl  4—  I 

S4-1 
while  1 

I  el  j  * —  0 

cl  j«-0-MPa 


e1 

e2 


count 


count4-^  count' 


ex  P,al,E  s 
e  x  P,bl -|-3,e 
tan(P) 


^yJcount^U  Astarcount  ]i(|] 


El  count 

cI  count4- 

A  starC0Unt 

^count 

es 

Jl 

I  al*— al  -(-2  if  slcoum»ex  p,al,E 

bl*-bl  +  l  if  elcounpe  x  p,bl-j-3,E  s 

count*-  count  +■  1 
break  if  al -f~bl=6 
break  if  al>3 
break  if  bl  >  3 

tempsto re*- augment  e! 


36*1 

e  rubber  ’  j 

iii  :=1 ..  rows  e  mbber 

a  rubberjji 

Q  rubber’ 

E  rubberjji 

0 

0  J. 

plot ‘-stack  plotl, augment  e  mbber> 


°  rubber 
MPa 


Tensile  Flow  Stress  Determination  cu^x)  :=iimerp  Piot<1>,Piot<2> ,x 


,;= _ ,__f 

plot  1  rows(plot) 


plot  r0Ws(  plot) 


curve(x)dx-MPa 


a  t  =  53.354°MPa 
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Calculation  of  Degraded  Laminate  Stiffness  (Compressive): 

The  first  ply  failure  strain  is  first  calculated.  Subsequently,  a  progressive  failure  of  plys  is  allowed  and  the 
stiffness  calculated  based  on  the  ply  damage  history. 

(5  !=180deg 
plotl  !=  al4-I 
bl4-l 
countf-2 
nl«-l 

54—  1 

while  1 


al  j4-0*MPa 

acount,  l*“al 
acount,2<_bl 


I  El  min 
count 


€  x  P  ,  a  1 ,  £  <j 
E  X  P*^l  “H^»E  s 


|  ^count4-  E 1  count 'tan^  P  ^ 


‘  S‘arcojm  M  M  h  mbber 
N  0‘S+N  90 *s H — ; — r - 


fVI  _  ,  ^  l1  rubber  - 

**  N0*  Qnewajj  S+N90*  Qnew^j  ^  s"l - ^ —  Q  rubber 


n  y  1  count*”"  [L  starcount  J  j 
el  count 

^  count*”  ^  s^arcount  ^count 


al4-al-|-2  if  Elcounp  ex  P.al,ss 
bl4-bl-|-l  if  Elcount=  E  x  0,bl -|-3,e  s 

count4~count-|- 1 
break  if  al-f-bl5^ 
break  if  al>3 
break  if  bl > 3 

tempstoree- augment  e],_!1L 


£  rubber  • 


r  06i 

z  rubber^ 

,  iii  *~1  ■■  rows  e  jobber  v  rubbery  *” 

Q  rubber' 

0 

L  J 

0  J 

plot! -stack  plotl, augment  e  rubber*- 


Compressive  Flow  Stress  Determination 


curve(x)  t=Hnterp  p]ot<*>,plot<‘>,x 


P,ot<1>  rows(plot) 


plQt  *  rows(  plot)  •  0 


curve(  x)d  x-MPa 
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Force  vs  Axial  Crush  Distance  Determination 


nn  !-30  bbb  *=1 ..  nn 


1 .  Establishing  Displacement  Array  ,_Dmodei<bbb-  >) 

6bbb  '  - - 

2.  Calculate  Open  Section  Correction  Factor 

®(8)  •  curvature —  b 


4>( 6 )  i-2-asin 


J- 


r  o  “t*  d  cutout^0^  to )  —  a(  5 ) 
cos(  Ct)  j 


-  d 


cutout- 


1 (W 


r  o“l“^  cutout^  )“  #(  S ) 
~  COS(CD  ) 


/cw_.  n  cutouts^5) 

E(5)-1 - 2rn - 

3.  Calculate  Unbalance  Flow  Stress  Correction  Factor 

2-o  t 

a  := - _  a  =1.728 

°ct°t 

4.  Calculate  Crushing  Force  for  Displacement  Array 


1  3  1 

■ 

i  3  r 

1 

3  r 

pbbb:=if 

_  77  *7 

a  Sbbb  c*  1  L  ’  5bbb 

a  8bbb  8 bbb  ' 

4  n*j2-c  c-a2  tLT-  8bbb  2 

,4-7r aJi>g  c-a^ -I 

1  c  1 

[  L  •  6bbb 

PI(xx)  :=( linterp( 6 ,P,xx)) 

Energy  Absorption  of  the  Structure: 


Energv(61 1 )  ’ — linterp( 5 ,P,61 1 ) 


Cumulative_Energ>j)kk 


rbbb 

m 


Energv(  x  )d  x 


Cumulative_Energyrows^  g  ^  =  0.058°MJ 
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Appendix  B 


Frustum  Section  Analysis 
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Ply  Stiffness  and  Layup  Inputs: 


1 .  Material  System:  E-Glass  (1583)/Epoxy  (7780)  woven  cloth 

2.  Layup  Configuration:  [90, 45,0, -45k, 90, 45, 0,0, 45, 90, [90, 45, 0,-45]2 


N  o’=6  N  n45  1=4 

N9o:=6  N  p45  :=6 

3.  Total  Number  of  Plys  Total _piys:=  n  0+n  n45-i-N  p45+N  90 

4.  Ply  Thickness  h:=oi8in 

5.  Laminate  Thickness  t  L  :=Totai_piys-h 

L  t  l =0.0101*m 

6.  Strain  Applied  to  the  Material  (circumferential): 

p=atan(Sy/£x)  p  :=89.5<ieg 

Shear  Strain:  es:=o 

7.  Number  and  Orientation  of  Ply  Angles: 


Plys:  =4 

“  0  ■ 

a  “1 ..  Plys 

c  :=1  ..3  Plys 

e  := 

90 

45 

.~45. 

<Jeg 


8.  Engineering  Constants: 


ex:=3-msi 
E  y  :=2.7  MSI 


Es:=.S5MSl  Gm:=1.26GPa 

E  m  :=3.4  GPa 


sequence 


“90 

45 

0 

-45 

90 

45 

0 

-45 


90 

45 

0 

0 

45 

90 


•deg 


-45 

0 

45 

90 


-45 


0 

45 

90 


v  x  *=-15 
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Component  Parameters: 

1 .  Radius  of  Shell: 

R  shell :  =3.25  ft 

2.  Length  of  Shell: 

L  shell 

3.  Wave  Length  of  Elastic  Buckling  Wave: 

H  :=V|/R  shell1  L 

4.  Number  of  Buckling  Waves  Along  Shell 

,_L  she|i 

Ply  Strength  Parameter  Inputs: 

Number  waves  -3  H 

1.  Tensile  Strength  in  Longitudinal  Direction 

X:=39.1.KSI 

2.  Compressive  Strength  in  Longitudinal  Direction  v  .  „CI 

A  prime* 

3.  Tensile  Strength  in  Transverse  Direction 

y:=39.iksi 

4.  Compressive  Strength  in  the  Transverse  Direction 

5.  Shear  Strength 

Y  prime  :=33  8KS1 

S  :=9.4I  KSI 


Matrix  and  Fiber  Damage  Inputs 


^  f_star  *  -01 

n  y:=.5i6i 

n  !=.l 

E  m_star ’“-l5 

n  s  := 316 

F  xv_star 

Damage  Ply  Calculations: 


-  -Ex 
Efx-- 


Gfx:= 


,  v  m 

v  ml 

-  1 

,  vm 

v  m 

1-M  c' - 

1  5  Vf 

e  a.  := 

'+hyV7 

Es 

<5m  J 

E  IV 

ty 

^  m 

1-1 


E  v  damaserf  ^ )  * 


1 

v  m 

1  ,  1  n  y~ 

i+ny_ 

E  fy  F  E  m 

r  1 


E  s_damaged(  F )  • 


1 

v  m 

1  d  1  n  S  ’“ 

1  .  v  m 

1  H-n  s. - 

*  V  f 

r  1 
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Building  Ply  Data: 

Undamaged  ply,  ply  with  matrix  damage  and  ply  with  fiber  damage. 


E  xx  •=! 


'  Ex  ' 

Ey 

E  s 

E  X 

Eyy:  = 

E  y_damaged  E  m_star 

II 

*k 

W 

E  s_damaged  E  m_star 

E  x  'E  f_star^ 

E  y_damaged  E  f_star 

E  s_damaged  E  f_star 

vx 

E  xy_star 

vxx:= 

v  x'E  m_star 

EE  xy_star*~ 

E  xy_star'E  m_star 

v  x'E  f_star  _ 

E  xy_star'E  f_star 

xxq:=x 


^  primeq  prime' 


YYq:=Y 


^  primeq  prime 


ssq  :=s 


vyyq*=v  xx( 


yyq 

xxn 


Rotation  Matrices: 

ml 

cos 

n ' 

:=sin  0„ 

a 

a 

a 

a 

ml  2 

n  2 

a 

a 

a  a 

j  := 

n  2 

ml  2 

-2  ml  n 

a 

a 

a 

a  a 

-ml  -n 
a  a 

ml  -n 
a  a 

ml  2-  n 
a  a 

Rotation  la  H 


ml  4 

n  4 

2-  ml,  2- 

n 

2 

4 

ml 

2  n  2 

a 

a 

a 

a 

a 

a 

n  4 

ml  4 

2  ml,  2- 

n 

2 

4 

ml 

2-  n  2 

a 

a 

a 

a 

a 

a 

ml  2-  n  2 
a  a 

mla  2-  "a2 

"V+ 

na 

4 

-4 

mla 

2’  na  2 

ml  2*  n  2 

ml  2-  n  2 

-2  ml,  2 

•  tu 

2 

f  ml  2- 

-  »  2f 

a  a 

a  a 

a 

a 

L 

a 

3  J 

ml  ^-n 

3 

-ml  •  n 

ml  ■  n  ^  — 

ml 

3-n 

:i 

ml  • 

n  3 

-  ml. 

a  a 

a  a 

a  a 

a 

l  a 

a 

a 

a 

ml  ■  n  ^ 

-  ml 

ml  ^ -n  —ml  • 

n  3 

2\ 

ml 

3n  ■ 

—  ml  •  n 

a  a 

a  a 

a  a 

a 

a 

a 

a 

a  a 

Rotation2a  !=j 


ml 


ml 


mlana  -m,a"a 
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Ply  Stiffness  Calculation: 

Calculates  the  zero  direction  stiffness  matrix  for  the  various  ply  states  (undamaged,  matrix  damaged,  and 
fiber  damaged) 


l-v  . 


V»q 

*q'Eyyq 


l-v 


xxq  yyq 


o 


vyyqExxq 

|-vxxqvyyq 


E 


yyq 


0 


0 


Creates  the  stiffness  matrices  for  every  combination  of  ply  angle  and  ply  state.  Qnew  is  a  3  x  4  matrix 


Q 


r°tq,a 


r  Q 

1,1 

Q 

vq  2,2 

^rolq,a  i  ^rolq,a  3  ^rotq,a  5 

-Rotation  1  • 

- 

a 

Qq  1,2 

Qnewq,a* 

^rolq,a  3  ^rotq,a  1  ^rotq,a  g 

Qq  3,3_ 

Q  r°tq , a  5  Qrotq,a6  Qrotq,a4 

Flexural  Stiffness: 

1.  Indexing  variables  fr:=i..Totai_piys+i 

gg  !=1  ..Total_plvs 

2.  Establishing  an  array  of  thickness  positions  for  each  ply 


Zfl.-llt+tff-Dh 


3.  Assigns  the  appropriate  stiffness  matrix  based  on  the  layup  sequence 


Q  flex 


gg 


Q  new)  >  |  sequencegg=0<ieg 
Q  newj  ^  if  sequencegg=90deg 
Qnew13  'f  sequencegg=45-deg 
Q  newj  4  otherwise 
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Stress  Space  Strength  Parameters: 

A  separate  F  must  be  calculated  for  each  ply  state 


f*  := 


FF 


vv  ,vv'  ", - 

AAq  AA  pnmeq 


xy_starq 


FF 


xy_starq 


4 


XXqXXprimei).YYq.YYprimeq 


^XXqXX  pnmeq 'YYqYY  prim£q 


YYq'YY  primeq 


SS„ 


~  1  1 

_  primeq 

1  1 

Y  q  Y  Y  Primeq 


Strain  Space  Strength  Parameters 

• 

• 

1 .  Strain  strength  parameter  for  each  ply  state 

.  .0  0  ... 

^1,1  a  1,1  q  1,2 

0 

+2  Faq  1,2  Q<i  l.i  Qq  1,2  ''  + 

Faq  1,2 [  Q<i  i.r  2,2+[  Qq  1,2]  ] 

... 

+  Faq  2,2'[  Qi  1,2]  + 

F,  •  Q  •  Q, 

*<1  2,2  4  1,2  4  2,2 

Gbaselq:“ 

Faq  i.r  Q(i  i.r  Q<i  1,2"' 

0 

+  Faq  1,2[  Qq  u  'Q<i2,2+[Qq  1.2 f\~ 

+  2*  Faq  1.3  0l)  1.2  2.2*** 

+  F  0  •  Q  •  Q 

®q  2,2  q  i ,2  q  2,2 

0 

0 

FX,  3.3  [ 

G  base2q  ;=| 


%,•  Qq 

F  bq  | '  Qq 


i,l  +  F bq  2  ^  1,2 
1,2+  F  bq  2  2,2 
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2.  Calculating  G  matrices  for  all  combinations  of  ply  state  and  angle: 


G  l  trans- 


Gbase2j  ^ 

Rotation2, 

i 

Gbase2]  2 

G  base22  j 

Rotation2, 

G  basel]  j  j 

1 

G  base22  2 

G  basel]  ^ 

Gbase23  ^ 

Rotation  1. 

Rotation2, 

1 

G  base  1  ]  j  2 

1 

Gbase23  2 

G  basel j  ^  ^ 

Gbase2j  ^ 

G  basel 2  j  ] 

Rotation^ 

Gbase2j  2 

G  basel2  2  2 

G  base22  j 

Rotation  lj 

Gbasel2  j  2 

Rotation22 

G  base22  2 

Gbase!2  3  3 

Gbase23 

G  basel3  ^  ^ 

Rotation22 

Gbase23  2 

G  basel 3  2  2 

G  2trans 

Gbase2j  ^ 

Rotation  1  j 

G  basel 3  j  2 

Rotation^  • 

Gbase2j  2 

G  basel 3  33 

G  base22  j 

G  basel j  j  j 

Rotation^  • 

G  base22  2 

G  basel j  2  2 

Gbase23 

Rotation^* 

G  baselj  j  2 

Rotation^  • 

Gbase23  2 

G  basel]  ^  ^ 

Gbase2j  ^ 

I  ’  2 

Rotation2r 

G  basel 2  j  j 

4 

G  base2j  2 

G  basel  2  2  9 

Gbase22  j 

Rotation  ^  • 

Rotation2.  ■ 

G  basel2  j  2 

4 

G  base22  2 

G  basel?  3  3 

Gbase23  ^ 

;  ’  * 

Rotation2,  • 

G  baselj 

4 

G  base23  2 

Gbasel3  2  2 

Rotation  lv 

X 

Gbasel3  (  ^ 

G  basel 3  ^  3 

G  basel j  ,  , 

G  baselj  -  ^ 

Rotation  U  • 

X  ,  X 

j 

G  basel]  ,  _ 

j  j  2 

G  basel |  3  3 
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Rotation  U 

L 

Gbasel2  , 

G  basel2  2  2 

j 

Rotation  U 

Gbasel2  j  2 
Gbasel2  ^  ^ 
Qbasel3  j 

G  base! 3  2  2 

j 

G  basel3  2 

Rotation!.  • 

G  base  1 3  ^  ^ 

'  GbaselI  u' 

G  basel |  2  2 

4 

G  base| i  I>2 

G  basel |  ^  ^ 

Rotation  I.  * 

Gbasel2  j  ^ 

G  basel2  2  2 

°basel2  ,  2 
Gbasel2  ^  3 

Gbasel3  j  ^ 

Rotation  L- 

Gbasel3  2  ^ 

4 

G  basel 3  ^  2 
Gbasel3  3  3 

G  ltransj.  j  G  ltrans^  3 
G  ltran^  ^  G  Itran^.  n 
G  ltrans^.  .  G  ltran^,  6 


G  ltran^  ,  G2trans^  j 

G  ltrans^  6  G2tran^  2 

G  ltran^.  4  G2trans^  ^ 


Transverse  Strain 


Formation  of  Strain  Space  Failure  Envelopes 


1 .  The  following  relationship  is  used  to  determine  the  failure  envelope:  Gi .jcrEj'l-Gi'ci“1 

2.  After  shifting  to  polar  coordinates,  the  above  equation  can  be  solved  for  r  using  the  quadratic  equation: 

S|*r-cos$  £-,=r'sin<t>  A-i^-f-B-r— 1*0 

3.  Calculation  of  coefficients  for  the  quadratic  formula 

A(<t>  ,b)  :=  Gb  2-sin(ij))2-t-2-  Gb  ]  ^-cos< 4> )-sin( <f>  )-|-  Gb  (  t  -cos( <j>  )2 
B  <(> ,b,e  s  :=2-  Gb  ■sin(+)-es+.  Gb  4-sin(<|>)+2-  Gb  (  3«os(<H-e  s+  Gb  ,  4<»*(*) 

C  <t> . b,e  s  •”  Gb  3  4‘£  s+  Gb  3  3  £  s2—  1 

-B  *,b,Es  +Jb  ,b,e  s  2-4  A(*,b)  C  (fr.b.E,; 

r  i,b,£c  !- - 2 - 

2  A(  <> ,  b ) 

Strain  Space  Failure  Envelopes: 

<(>  :=0<Jeg,  1 ‘deg..  360  deg  y(x)  :=tan(p)-x  x  :=0,^2f^2..cos(p)  £x4>,b,Es  t=r  <t>  ,b,e  s  -cos{  4> )  Ey$,b,Es  l=r  ^,b,E  s  sin($) 


Ply  Strain  Strength  Envelopes 


Longitudinal  Strain 

-  0  ply-Undamaged 

.  0  ply-matrix  damaged 

.  90  plv-Undamaged 

-  90  plv-matrix  damaged 

-  45  ply-Undamaged 

45  ply-matrix  damaged 

- .  -45  ply-Undamaged 

-  -45  ply-matrix  damaged 

-  Strain  Relationship  Line 
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Calculation  of  Degraded  Laminate  Stiffness: 

The  first  ply  failure  strain  is  first  calculated.  Subsequently,  a  progressive  failure  of  plys  is  allowed  and  the 
stiffness  calculated  based  on  the  ply  damage  history. 


plot 


al«-l 


bl*-1 

cl«-I 


dl  4—  1 


count*- 2 
nl«-l 
while  1 

Icl  |  *—  0 

ol  j4-0  MPa 


el 


count 


min 


^count4”8*  count 


e  x  (3,bl h-3,e  s 

E  x  P,Cl-f>6,£  s 

e  x  p,dl  9,e  s 
■tan(P) 

1 


A  i 

starcount^l^o^N 90+N  n45+N  p45 

.-1 


^  y^count4-^  ^  slarcount  J2 , 2J 

A  starcount 


Qnewal  ,  90*  Qnewbl>2  +Nn45*  Qnewcl|3  +N  p45*  Q  newdl  ,4  ] 


a! 


count 


’si 


count 


e2, 


count 
£  s 


al«-al  +  l  if  £lcount=Ex  MUS 
bl<-bl+2  if  El^unpEx  P,bi-|-3,e  s 

cl<-cl-H  if  ^count^x  P.cl+6,es 
dl<-dl  +  l  if  e1count*ex  p,dl+9,es 

count*- count-}- 1 

break  if  al +bl -}-cl-f-dl=12 

break  if  al>3 

break  if  bl>3 

break  if  cl>3 

break  if  dl>3 

*  ol 

tempstore*- augment  e2,-— - 


augment  tempstore 


H  yl 

’TjET 
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Ply  Transverse  Young's  Moduli  as  a  function  of  Strain 


Transverse  Strain 


0 

0 

0 

0.014 

247.801 

16.449 

0.015 

180.252 

10.741 

0.016 

126.253 

7.07 

0.025 

85.233 

2.765 

0.171 

94.954 

0.554 

Transverse  Stiffness  Variations  with  Axial  Location  and  Axial  Crush 
Distance: 


1.  Setting  up  stations  along  the  x'  axis  and  axial  crush  distances  at  which  Ey  is  to  be  calculated 

n  !=2000  yyy‘=I..  17  zzz ’=1 .. n 


k  Pr,mezzz 


'=  222  1  u  „ 

•  “n  Tii  ,Heff 


1  temp  *' 


2-10* 


r3 


r  3 


3-10* 

510“ 


610 


r3 


w-3 


7-10* 

8  •  1 0“  3 

910“3 

0.011 

0.014 

0.02 

0.029 

0.041 

0.068 

0.093 

0.1 

.2 


Heff 
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2.  Calculation  of  circumferential  strain  based  on  axial  position  and  crush  distance 


._x  prim e^-5  tempyyy  j  Heff  “7 

Ee^yyy  '  r  shell*  ar  *J°  temp^,"7 

3.  Assignment  of  the  appropriate  Young's  moduli  to  the  correct  axial  position  and  crush  distance  based  on  the 
transverse  Young's  Moduli  vs  strain  plot  and  circumferential  strain. 


OGPa  if  e  a  >171 

°zzz,vyy 

.554-GPa  if  .171>e  *  >025 

°zzz,yyy 

“0 

0 

0 

0.014 

247.801 

16.449 

2.765-GPa  if  .025>s  >016 

°zzz,yyy 

0.015 

180.252 

plot  = 

10.741 

7.07GPa  if  .016>e  fl  >015 

°zzz,yyy 

0.016 

126.253 

7.07 

10.741GPa  if  .O15>e0  >014 

°zzz,yyy 

0.025 

85.233 

2.765 

0.171 

94.954 

0.554 

1 6.449  GPa  otherwise 

Flexural  Strain  and  Stiffness  Calculations,  and  Thickness  Variations: 

1.  Flexural  Strain  as  a  function  of  crush  distance  .oi8  in  =  4.572»io-4*m 


e  xf<6.gg) :_2' 


Wl+2' 


2R  shell 


R  shell1  L  — 4- 


8  :=0,.00I..2v^R  shell  t  L 


Plv  5 
Ply  9 

Tensile  Failure  Strain  for  0  Oriented  Ply 


Flexural  Strain  F|exura,  Strain 


2.  Determination  of  flexural  stiffness  and  thickness  variations  with  crush  distance.  This  a  manual  operation  where 
the  plys  are  failed  progressively  as  crush  distance  is  increased  to  2H  using  the  graphs  above  and  the  included  ply 
strain  limits.  The  thickness  is  determined  by  removing  ply  thicknesses  as  they  fail. 


'  star  * 


12 


(2-h) 


I 

rr 


Total_pIys 

Z 1  Q  flexxx([z(  xx+ 1  )]3_  \x  3]  ssxx 


-1 


E  xf  *  d  starj  j 


E  xf  =  20.6843°GPa 


■  0 

.002 

'16.4811* 

"22h’ 

.003 

16.4811 

22  -h 

.005 

17.3293 

19-h 

.006 

15.9983 

15  h 

.007 

16.4377 

13  h 

.008 

18.3322 

1 2-h 

.009 

16.5229 

1 1  -h 

.011 

H  efT 

16.7947 

10  h 

.014 

'  ft  ‘m  ^flexure*” 

17.0984 

•GPa  t  f  := 

9h 

.02 

14.2239 

8h 

.029 

16.142 

7-h 

.041 

17.2425 

6h 

.068 

18.3592 

5-h 

.093 

12.799 

4h 

H 

18.4517 

3h 

rn 

20.6843 

2h 

2-H 

20.6843 

2-h  m 

m 
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Force  vs  Axial  Crush  Distance  Determination: 


1 L  2  c 

pyyy '-2’7t  R  „  2'  2-  Eeppp,yyy’  xPnmeppp 
K  shell"  eff  ppp=  I 


2  H  efT  E  flexure^*  shell'  ‘  fyvv  8 

_  h - : - - - —  •  Heff - -l 

"  3-Heff3  2 


Enerqv  Absorption  of  the  Structure: 


Energy(8U)  ;=linterp(8,P,811) 


|*8yyy 

Cumulative_Energ>^  1=  J  Energy(  x)d  x  Cumulative J:nergyrQW^  5  y  -  0.021  °MJ 
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Appendix  C 


Cylindrical  Section  Analysis 
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Ply  Stiffness  and  Lay-up  Inputs: 

1.  Material  System:  E-Glass  (1583)/Epoxy  (7780)  woven  cloth 

2.  Lay-up  Configuration:  [90,45,0,-45]3,90,45,45,90,[90,45,0>-45]3 


N0:=6  N  „45  :=6 

N  90 :=8  N  p45 :=8 

3.  Total  Number  of  Plys  Totai_piyS:=  n  0+n  n45+N  p45+N  90 

4.  Ply  Thickness  h:=.oi8in 

5.  Laminate  Thickness  t  L  :=Totai_piys-h 

6.  Strain  Applied  to  the  Material  (circumferential): 

P=atan(sy/£x)  p  :=89.5-deg 

Shear  Strain:  es:=o 

7.  Number  and  Orientation  of  Ply  Angles: 


Plys  :=4 
a  i=l ..  Plys 


c  :=1 ..  3-Plys 


8.  Engineering  Constants: 


HX:=3MSI 
Ey  :=2.7  MSI 


e:= 


*  0 
90 
45 
-45 


•deg 


Es:=.55MSI  Gm:=1.26GPa 

E  m  3.4'GPa 


‘90 

45 

0 

-45 


sequence 


90 

45 

0 

-45 

90 

45 

0 

-45 

90 

45 

45 

90 

-45 

0 

45 

90 


-45 

0 

45 

90 

-45 

0 

45 

90 


v  f  :=5 


Vx:=.i5 
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Component  Parameters: 


1.  Radius  of  Shell: 

2.  Length  of  Shell: 

3.  Wave  Length  of  Elastic  Buckling  Wave: 

(shell  only,  no  stiffener  effect) 

4.  Stiffener  Parameters 

Undamaged  Stiffness 

Number 

Area  of  Single  Stiffener 

5.  Open  Section  Parameters 

Reduction  Parameter 


R  shell  ■  3,5-ft 
L  shell  :=33  in 

H  :=^/R  shell1  L 

E  carbon  •~53.793GPa 
N  stiffener-** 

A  stiffener  :=(  5. 526  in  -f- ( 8  -in  -  2-.45in ) )  -.45-in 
C,  :=62.5-deg 

y:=iLz£  y  =  0.653 


Piy  Strength  Parameter  Inputs: 

1.  Tensile  Strength  in  Longitudinal  Direction 

2.  Compressive  Strength  in  Longitudinal  Direction 

3.  Tensile  Strength  in  Transverse  Direction 

4.  Compressive  Strength  in  the  Transverse  Direction 

5.  Shear  Strength 


x:=39.i-ksi 

X  prime -33.8-KSI 

Y:=39.1-KSI 


^  prime1  33.8-KSI 
S  :=9.41-KSI 


Matrix  and  Fiber  Damage  Inputs 


E  f_star 

n  y  := 5161 

n  *=.t 

^  m_star 

1  s  :=  316 

F  xy_star,—" 
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Damage  Ply  Calculations 


E»_  A 

fx 1  1 

V  f 


G  fx  H 


•  .  v  m 

v  m" 

11  S'T7 

-I  r 

V  I 

v  t 

Efy:= 

E  S 

G  m 

I+T1- 


E  y_damaged(F) 


J-hn  y 


v  m  [  fy 
vf 


Tly- 


+T 


E  s  damaged^  F )  ■ 


1-j-Tl  , 


*fe+T‘ 


vf 


Building  Ply  Data: 

Undamaged  ply,  ply  with  matrix  damage  and  ply  with  fiber  damage. 


'  Ex  ‘ 

Ey 

E  s 

^  XX  •“ 

E  X 

E  XE  f_star_ 

Ey>,:= 

E  y_damaged  ^  m_star 

E  y_damaged  ^  f_star 

Ess:= 

E  s_damaged  E  m_star 
b  ^  s_damaged  ^  f_star 

vx 

F  xy_star 

vxx:= 

v  x*E  m_star 

FF  xv_star= 

^  xy_star'E  m_star 

v  x  E  f_star  _ 

. F  xy_star‘E  f_star  _ 

xxq:=x 


xx 


primeq  prime'  -g- 


YYq.-Y  YYprimeq:=Y  prime 


Rotation  Matrices: 


mla:=cos0a  V'si"ea 


Ja:=l 


mla  2  "a  2 


2mlana 


na  2  mla  2  '2-mlana 


-ml  n  ml  -nn  ml  —  n 
a  a  a  a  a  a 


2 
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- 1 

3 

tu 

"a4 

2-  ml  2*  n  2 
a  .  a 

4  mla2-  na 

"a4 

mla4 

2  ml  2.  n  2 
a  a 

4  ml  2*  n 
a  a 

mla  2'  na 

2  mla2'na2 

ml  n  4 

a  1  a 

-4  ml,  2-  na 

Rotation  1  1- 

a 

mla  2-  na 

2  "I  2-„  2 
a  a 

-2  ml  2-  n  2 
a  a 

[  ml  2—  n  4 
[_  a  a 

m,a  3'na 

-m,a'  na  3 

mla-  na  3-  m1a  3-na 

-i 

mIa’  na  3~  ml; 

.  mla  na  3 

-  mla  3  "a 

3  3 

ml  a  -ml  -  n 
a  a  a  a 

2i 

ml  -ml  - 

a  a  a 

Rotation  != 

u 


ml  2  n  2 
a  a 

n2  ml  2 


ml  -n  -ml  n0 
a  a  a  a 


Ply  Stiffness  Calculation: 

Calculates  the  zero  direction  stiffness  matrix  for  the  various  ply  states  (undamaged,  matrix  damaged,  and 
fiber  damaged) 


<v=l 


vy>qExxq 


|-vxxq'vyyq  l_vxxqvyyq 
vxxqEyyq  Eyyq 


1  “  v  xxq  *v  yyq  1  ~  v  xxq  v  yyq 


Creates  the  stiffness  matrices  for  every  combination  of  ply  angle  and  ply  state.  Qnew  is  a  3  x  4  matrix 


Q  rnt  1  -Rotation  1 
r0lq,a  a 


Qn 

<n,i 

Q<i  2,2 

Qro‘q,a,  Qrotq,a3  Qrotq,a5" 

Q<1  1,2 

Qnewq,a: 

Qrotqaj  Qrotq  a  2  Qrotq  a  6 

Qq  3'3. 

^  rotq,a  5  ^  rotq,a  6  ^  rotq,a  4 

Undamaged  Elastic  Stiffness  of  Laminate: 

Q  laminate  -  £  ft  n45  +  N  ^  [N  O'  Q  new,  _ ,  +N  90' Q  new,  2  +N  n45'  Q  new,  3  -bN  p45‘  Q  new,  4  ] 


Q  laminate = 


17.575  5.256  -0.038' 
5.256  17.726  -0.038 

-0.038  -0.038  6.198 


°GPa  S  .-Q  |aminate 


p  —  1 
t  Eglass 

*1 
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Effective  Thickness  Calculations 


A  stiflener'N  stiffener  E  carbon 
(i’7t-2  <;)-R  Sheii  E  Eglass 


H  eff:=VR  shell 1  efT 


tefr=0.035«m  tL  =  0.0128*m 

H  eff=  0.194’m  H  =  0.117*m 


Number 


,_L  shell 
waves  ’  yjj 


eff 


Number  vvaves  =  2.159 


_  _  L  shell 

Number  waves  .=2  H  eff=0.21*m 

Flexural  Stiffness: 

1.  Indexing  variables  fr  :=i ..  Totai_piys-j- 1 

gg  *=1  ..Total_plys 

2.  Establishing  an  array  of  thickness  positions  for  each  ply 

3.  Assigns  the  appropriate  stiffness  matrix  based  on  the  lay-up  sequence 


Q  flex, 


gg 


Q  new,  1 1  'f  sequencegg=0-deg 

Q  newj  2  if  sequencegg=90-deg 
Qnew]  3  if  sequence  «45 -deg 
Q  new]  4  otherwise 


Stress  Space  Strength  Parameters: 

A  separate  F  must  be  calculated  for  each  ply  state 


1  FF  xy_starq 

^'^Prime<l  ^XXq.XX  pnmeq  YYq.YYpH^ 

xy_starq  I 

^XXq.XXprimeq.YYq-YYprimeq  YVq.YYprimeq 


’  1  1 
w-“  xx — • — 

AAq  AA  primeq 
1  _  1 
YYq  YY  primeq 


/ 


90 


Strain  Space  Strength  Parameters 

1 .  Strain  strength  parameter  for  each  ply  state 


G  baselq  H 


+ 

+  F. 


,Faq  1,1  Q<i  1,1’  Q<i  1,2"' 

!*  F  o  *  Q  •  Q 

“q  1,2  q  1,1  4  1,2 

+  Fflq  1,2^  Q<i  i.r  Qq  2,2+[  Qq  1, 

SwfVJ 

+  F-  *  0  0 

^  2,2  ^  1,2  q  2,2 

'  Qq  i.i  Qq  i,2"' 

fv,..[V2]2'" 

, 2_[  Q<1  l,rQ0  2,2+[  Qq  1,2]  ] 

+  2  Fafl  1.2  ^  t,2  Q1  2.2'” 

•  Q  •  Q 

2  q  1,2  q  2,2 

0 

0 

Ss.,{V3]\ 


-  F 

^base2q*”  p 


i,i  +  Fbq  2  Q(i  1,2 
1,2+  F  bq  2  ^  2,2 


2.  Calculating  G  matrices  for  all  combinations  of  ply.  state  and  angle: 


Rotation2j 

G  base2] 

G  base2j  2 

r 

G  basel]  j  j 

Rotation2^ 

Gbase22  f 

G  base22  2 

Rotationlj 

G  basel j  2  2 
G  base  1 j  j  2 

Rotation^  • 

Gbase23  j 

G  base23 

G basel,  3  3 
G  basel2  |  j 

Rotation^  ■ 

"  Gbase2j  j 

G  base2]  2 

Rotation  lj  • 

^  basel 2  -j  2 
^  basel 2  j  1 

Rotation22  • 

'  G  base22  1 

G  base22  2 

Rotation  I  j 

G  basel 2  ^  . 

"  Gbasel3  , 

G  baselj  22 

G  2trans  • 
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Rotation^  • 

Rotation^  • 

Gbase23  j 

G  base23  2 
'  G  base2]  J 

G  1  trans  • 


Rotation  lj 


I 

Rotation  L 

G  base  1 3  ^  2 
G  basel 3  3  3 
G  basel]  j 

G  basel]  2  2 

2 

Rotation  1- 

G  base> i  1>2 

G  basel ]  3  3 

Gbase,2  U 

G  basel 2  2  n 

Rotations 

cbasel2  ]>2 
Gbasel2  3  3 

G  basel 3  j  j 

G  basel 3  ^  ? 

L 

Rotation  U- 

G  basel 3  ,  2 

G  basel3  3  3 

G  basel i  j 

G  basel  ]  n  2 

3 

Rotation  U  ■ 

G  basel]  j  2 

G  basel]  ,  - 
G  basel2  j  j 

Gbaseb  2 

j 

Gbasel2  j  2 
Gbasel2  3  3 

Rotation  U- 

G  basel 3  j  } 

G  basel 3  0  2 

3 

Gbasel3  j  0 

G  basel 3  3  3 

Gbaseh  ,1 

Rotation  1 .  • 

1  1,1 

G  basel]  2 

4 

G  basel]  j  2 

G  basel,  3  3 

G  basel2  j  j 

Rotation  14* 

G  basel 2  n  0 

G  basel t  j  2 

G  basel 2  3  3 

G  basel3 

Rotation  1^- 

G  basel 3  2  2 

Rotation^ 

Gba$e2|  2J 

Rotation^ 

’  G  base22  ' 
0  base22  2 

Rotation^  ■ 

G  base23  ( 

0  base23  2 

Rotation24  • 

Gbase2)  j 
Gbase2,  2 

Rotation2^- 

Gbase22  j 

0  base22 

Rotation2^- 

r  G  base23  ' 
0  base23  2 
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Rotation^  • 


Gbasel3  j  2 
Gbasel3  ^  ^ 


G  Itrans^  j 

G  Itransj.  ^ 

G  ltran%  _ 

G  2transc 

Gc:= 

G  ltran^  ^ 

G  Itran^  2 

G  ltran^  6 

G  2transc 

G  ltran^  . 

G  Itransg  ^ 

G  ItranSj.  ^ 

G  2transc 

Formation  of  Strain  Space  Failure  Envelopes: 

1 .  The  following  relationship  is  used  to  determine  the  failure  envelope:  Gj  j  -Cj  -Ej  +-Gj  -e«i 

2.  After  shifting  to  polar  coordinates,  the  above  equation  can  be  solved  for  r  using  the  quadratic  equation: 

E]®rcos<|(  E2SSrsin<t  Ar2  fBr-  1=0 

3.  Calculation  of  coefficients  for  the  quadratic  formula 

A(<t>  ,b)  :=  Gb  2  i-sin((|i)2+2-  Gb  (  ^ -cos( <f> )  sin( 4> )  +-  Gb  f  ^  cos( <(.  )2 
B  <l> , b,e  s  :=2-  Gb  ^  3-sin(<(.)  E  s+-  Gb  ^  4-sin( 4> Gb  f  3 cos( <(> )  e  s+  Gfe  (  4<os(<|>) 

C<f,b,Es  :=  Gb  3  4es+  Gb  3  3-es2-1 

_-B  $,b,Es  +  WB  *,b,Es  2-4A(*,b)-C  <(> ,b,e  s 
r  di.b.E  ,  — _ 1 _ _ _ 
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-  0  ply-Undamaged 

.  0  ply-matrix  damai 

- .  90  ply-Undamaged 

-  90  ply-matrix  dam; 

-  45  ply-Undamaged 

45  ply-matrix  damj 

. 45  ply-Undamagei 

-  -45  ply-matrix  dam 

-  Strain  Relationship 


Calculation  of  Degraded  Laminate  Stiffness: 

The  first  ply  failure  strain  is  first  calculated.  Subsequently,  a  progressive  failure  of  plys  is  allowed  and  the 
stiffness  calculated  based  on  the  ply  damage  history. 


plot! 


al«-l 
bl«-l 
cl«-l 
dl  <—  1 
counts  2 
nl  4—  L 
i  while  1 
cl  j<-0 

ohf-OMPa 


l£lcount<-min 


£  X  P»a^ ’c  s 

e  x  p,bl+3,c  s 
e  x  P.cl+6,es 
e  x  p.dl-)-9,E  s 


fe2count<-E,coumtan(l}> 


*■  star. 


rcount<_  No-t-N9o+-Nn45+Np45  [N  °'  Qnewal,l  +N90'  °newbl,2  +Nn45'  Qnewc|3  +Np45'  Qnewdl  4  ] 
Ey1count<_[[  A  starcoUm  "'jj  2] 


cl 


count* 


slarcount 


‘cl. 


al*-al  +  l  if  elcoun=Ex  P,al,Es 
|bl«-bl-)-2  if  e1C0UM=e  x  p,bl-t-3,E  s 
cU-cl-H  if  elcount=£  x  P.d-t-6,Es 
dl«-dl-H  if  ''count®6*  P,dl+9,ES 
counts-  count-f- 1 
break  if  al -f-bl -}-cl-f-dl=12 
break  if  al  > 3 
break  if  bl > 3 
break  if  cl>3 
break  if  dl>3 

cl 


;empstore«- augment  e2, 


MPa 


luement  tempstore, 


Ev 


yt 

"CFa 
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Transverse  Stiffness  Variations  with  Axial  Location  and  Axial  Crush 
Distance: 


1 .  Setting  up  stations  along  the  x'  axis  and  axial  crush  distances  at  which  Ey  is  to  be  calculated 

yyy  !=1 ..  18  n  1=2000  zzz  !=1 ..  n 


zzz  1  u 

xpnmeZ2Z*  ~T~Tn  eff 


$  temp 


“  0  ' 
.002 
.003 
.004 
.005 
.006 
.007 
.01 
.013 
.015 


HefT 

_ m 


.02 

.027 

.028 

.043 

.06 

.115 

.155 

234m 


/ 


2.  Calculation  of  circumferential  strain  based  on  axial  position  and  crush  distance 


,=XPrime2Zz'6  temPy>'y  I  Hefr  ~ 

*  ezzz.yjy  Rshell  H  eff  '  J i>  temPyyy  “  4 

3.  Assignment  of  the  appropriate  Young's  moduli  to  the  correct  axial  position  and  crush  distance  based  on  the 
transverse  Young's  Moduli  vs  strain  plot  and  circumferential  strain. 


0-GPa  if  e  q  >171 

°zzz,yyy 

.469  GPa  if  ,I71>e  6  >025 

uzzz,vyy 

“0 

0 

0 

0.014 

245.846 

16.154 

2.779-GPa  if  025>s  >016 

° zzz, yyy 

0.015 

plot  = 

174.15 

10.133 

6.094GPa  if  ,016>e  8  >015 

°zzz,yyy 

0.016 

113.553 

6.094 

10,l33GPa  if  ,O15>E0  >014 

0.025 

90.216 

2.779 

° zzz,  yvv 

0.171 

80.302 

0.469 

16. 154  GPa  otherwise 
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Flexura!  Strain  Flexural  Strain 


Flexural  Strain  and  Stiffness  Calculations,  and  Thickness  Variations: 


Crush  Distance 


-  Plv  2 

.  Plv  4 

-  Plv  6 

-  Ply  8 

—  Ply  10 

.  Tensile  Failure  Strain  for  45/-4S  Oriented  Ply 

— — -  Compressive  Failure  Strain  for  45/-45  Oriented  Ply 

-  trace  8 

-  trace  9 


2.  Determination  of  flexural  stiffness  and  thickness  variations 
with  crush  distance.  This  a  manual  operation  where  the  plys  are 
failed  progressively  as  crush  distance  is  increased  to  2H  using 
the  graphs  above  and  the  included  ply  strain  limits.  The 
thickness  is  determined  by  removing  ply  thicknesses  as  they  fail. 


TotaI_plys 

Dstar:=“3  l'  Qflexxx{ftxx+I)]3-  *xx  3]  ssxx 

d  star* =D  star  1 
E  xf 1  “  d  starj  j 
E  xf=  1 1.3874°GPa 


V 

.002 

.003 

’16.4023’ 

’28  h* 

16.4023 

28-h 

16.7798 

25 -h 

.004 

15.9182 

20 -h 

.005 

17.3134 

18h 

.006 

16.0961 

16h 

.007 

1 5  -h 

.01 

16.4632 

16.0174 

13-h 

.013 

Heir  c 

16.5986 

•GPa  t  f  *— 

1 1  h 

.015 

.02 

■  . E  n  " 

jq—  111  tiexure 

19.1109 

10  h 

16.703 

9h 

.027 

17.9956 

8h 

.028 

17.8391 

7h 

.043 

.06 

13.7589 

6-h 

!  16.6569 

5h 

.115 

17.9459 

4-h 

.155 

23.1851 

3-h 

2-H 

[l  1.3874_ 

2h 

m 

■  J 
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Appendix  D 
Stiffener  Analysis 
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Ply  Stiffness  and  Lay-up  Inputs: 


1.  Material  System:  Carbon,  woven  cloth  (6K  5HS  satin  weave  IM7) 

2.  Lay-up  Configuration:  [90,45,0,-45]4[90,45,0,-45]4 


N0:=8  N  n45  :=8 

N  90 :=8  N  p45 :=8 


3.  Total  Number  of  Plys  Totai_piys:=  n  0-(-n  n45-t-N  p45-i-N  90 

4.  Ply  Thickness  h  :=.oi3s-in 

5.  Laminate  Thickness  t  L  :=Totai_pivsh 

6.  Strain  Applied  to  the  Material  (circumferential): 

p=atan(sy/ex)  p  :=89.s-deg 

Shear  Strain:  t  s  :=o 

7.  Number  and  Orientation  of  Ply  Angles: 


Plys  :=4 
a  1=1 ..  Plys 
c  :=1 ..  3  -PIys 

8.  Engineering  Constants: 


•deg 


EX:=11.6MS1  Es:=.37MSI 


E  y  :=1 1.6MSI 


E  m  :=3.4  GPa 


1  m  •" 


vf:=5 


vx:=  124 


1.26  GPa 
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Component  Parameters: 

1.  Radius  of  Shell:  r  sheii  :=3.5  ft 

2.  Length  of  Shell:  .  ... ... . 

S  L  shell  -4  471ft 

3.  Wave  Length  of  Elastic  Buckling  Wave:  _ 

H  shell'1  L 

4.  Number  of  Buckling  Waves  Along  Shell 

Number  waves 


Ply  Strength  Parameter  Inputs: 

1.  Tensile  Strength  in  Longitudinal  Direction 

2.  Compressive  Strength  in  Longitudinal  Direction 

3.  Tensile  Strength  in  Transverse  Direction 

4.  Compressive  Strength  in  the  Transverse  Direction 

5.  Shear  Strength 


x:=39.iksi 

X  prime -33.8KSI 

y:=39  i-ksi 

Y  prime  :=33.8  KSI 

S  :=9.41KSI 


Matrix  and  Fiber  Damage  Inputs 

E  f_star  T|  y  .5161  n.— .1 


Damage  Ply  Calculations: 


E  fx  :_T 


*  fx  ’ 


'+ns  — 

v  m 

_ns'“ 

- 1 

E  a.  := 

,  v  m 

l+r,  y._ 

v  m 

*\  v- - 

V  V  f 

E  s 

^  m 

Ey 

E  m 

^  y_damagerf  E ) 


E  s_damaged(  E ) 


1 

'  +1, 

v  m  1 

n  y~ 

Jfy  F 

Em  J 

I 

_L+i> 

v  m  11 
x\  s- - 

5  vf 

,  v  m 

'Tls  — 

Gfx  +  F 

Gm  j 
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Building  Ply  Data: 

Undamaged  plv.  ply  with  m 


Undamaged  ply,  ply  with  matrix  damage  and  ply  with  fiber  damage. 


XX 


v  x  -E  m_star 
v  x  E  f  star 


FF 


xv  star' 


E  xy_star‘E  m_star 
[  E  xy_star'E  f_star 


'  Ex  • 

Ey 

Es 

E  xx  *“ 

Ex 

Eyy:= 

E  y_damaged  E  m_star 

Ess ,= 

E  s_damaged  E  m_star 

E  x  'E  f_star_ 

E  y_damaged  E  f_star 

E  s_damaged  E  f_star 

r 

vx  I 

f  E  x\_star  1 

xx( 

?:=x 

®  primeq  prime'  — — 


YYq:=Y 


YY  primeq  ^  prime 


ssq:=s 


,  -=v  »'q 


Rotation  Matrices: 


mla:=cos  6a 


"a :=sin  6 a 


V=l 


ml  2  n2 


2-mlana 


"a  ‘  m,a  -2  mlana 

-mlana  mla-na  mla  2-  na  2 


Rotation  1 

a 


ml  4 

n  4 

2 •  ml,  2- 

n  2 

4 

ml  2-  2 

a 

a 

a 

a 

a  a 

n  4 

ml  4 

2  mia  2- 

n  2 

4  . 

ml  2*  n  2 

a 

a 

a 

a 

a  a 

ml  2*  n 
a  a 

2 

*)  ? 
m,a  •  na 

mlaV 

"a4 

-4 

ml  2'  n  2 
a  a 

ml  2-  n0 

2 

ml  2-  n  2 

-2  ml  2 

•  n  2 

r  ml  2-  n  212 

a  a 

a  a 

a 

a 

L 

a  a  j 

ml  3n 

-ml  *  nQ  3 

ml  •  n  3- 

ml  3*na 

2-fmL- 

n  3—  ml  3-i 

a  a 

a  a 

a  a 

a  a 

L  a 

a  a 

ml  •  n0  J 

-  ml  3n 

3  3 

ml  -n  —  ml  •  n 

2-["  ml 

3 

•n  —  ml  •  n 

a  a 

a  a 

a  a 

a  a 

L  a 

a  a  a 

Rotation2a  !=j 


ml 


ml 


ml  'n,  -ml 
a  a  a  a 
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Ply  Stiffness  Calculation: 

Calculates  the  zero  direction  stiffness  matrix  for  the  various  ply  states  (undamaged,  matrix  damaged,  and 
fiber  damaged) 


<v=l 


|-vxxq'Vyyq 
v  xxq  Eyyq 


v  »'q 


v  y>'q:E  xXq 
|-vxxq'vyyq 

Eyyq 

1-v  xxq'v  yyq 
0 


0 


0 


Creates  the  stiffness  matrices  for  every  combination  of  ply  angle  and  ply  state.  Qnew  is  a  3  x  4  matrix 


Hq,a 


'Vi' 

2,2 

"  Qrotq,a  | 

Qrotq,a  3 

Qrotq,a  5 

=Rotationl  • 
a 

Qq  i,2 

Qnewq  a,= 

Qro,q.a  3 

Qrotq,a  2 

Qrotq,a  6 

.  Qq  3'3. 

^  rotq,a  5 

Qrotq,a  6 

Qrotq,a4_ 

Undamaged  Elastic  Stiffness  of  Laminate: 


Q  laminate  N  0+N  90+N  n45+N  ^ 

0'  Q  newj  ^  j 

64.715 

26.586 

6.322*1 0" 15  " 

Q  laminate** 

26.586 

64.715 

-6.322*10" 1S 

6.322*10"15 

-6.322*1  (T^ 5  19.064  1 

Q  newi  2  H"^n45‘ 


^  ’“Q  laminate 


Q  ne\vj  3  p45-  Q  neWj  ^  J 


53.793°GPa 


53.793°GPa 
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Appendix  E 


Forward  Hemi-head  Pressure  Hull 
Analysis 
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Analysis  of  HY-80  Forward  Pressure 
Hull  Hemihead 


Stress  Strain  Relationship  for  HY-80 


Material  Flow  Stress  Determination: 


,•8 


a  r  :=- 


"  rows(  e  )  4 


rows(  e ) 


curve(x)dx 


a  c  =517.499 oMPa 


Structural  Properties: 


ii . 

1 L  '=16  m 


a  :=  1 


bbb  :=  1 ..  rows(8) 


Crushing  Force  Prediction: 


3  2 

Pbbb:=4  ,t  V2-Oc-0t-tLI-  5bbb  2 

Energy  Absorption: 


Energy(Sll)  :=  linterp( 5 ,P,5l  1 ) 


8  := 


'0.028 

0.1 

0.241 

0.242 

0.245 

0.247 

0.248 

0.25 

0.252 

0.254 

0.259 

0.264 

0.268 

0.277 

0.289 

0.291 

0.318 

0.349 

0.447 

0.519 

0.66 


m 


Cumulative_Energybbb  := 


•5bbb 

0-m 


Energy(x)dx 
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